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ABSTRACT 


The  multistatic  radar  offers  many  advantages  over  monostatic  radar  in  certain 
applications,  especially  since  the  receiving  stations  may  be  located  at  covert  and  distant 
sites  relative  to  the  transmitting  stations.  Furthermore,  continuous  wave  radars  are 
relatively  simple  and  inexpensive  to  employ  and  maintain.  Hence,  the  impetus  for 
developing  a  CW  multistatic  radar  system  for  high-resolution  imaging  was  conceived. 

This  thesis  is  a  proof  of  concept  demonstration  that  a  Doppler-only  multistatic 
radar  system  can  be  employed  to  provide  high  resolution  imaging  of  airborne  targets  in 
support  of  non-cooperative  target  recognition.  Through  an  understanding  of  conventional 
imaging  techniques  and  formulation  of  the  inverse  problem  in  radar  imaging,  a 
demonstration  radar  model  based  on  one  transmitter  and  two  receivers  was  designed  to 
determine  the  accurate  position  and  velocity  of  simulated  targets.  The  extraction  errors 
resulted  from  the  range,  bearing  and  velocity  measurements  were  congruent  with  the 
physical  limitations  of  each  transmitter-receiver  pair.  Through  the  employment  of  a 
multistatic  system,  the  geometrical  diversity  allowed  these  limitations  to  be  overcome. 
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I.  INTRODUCTION 


A.  DOPPLER-ONLY  MULTISTATIC  RADAR  IMAGING 

IEEE  Standard  686-1997  [1]  defined  the  multistatie  radar  as  a  radar  system 
having  two  or  more  transmitting  or  reeeiving  antennas  with  all  antennas  separated  by 
large  distanees  when  eompared  to  the  antenna  sizes.  It  has  at  least  three  eomponents  — 
for  instanee,  one  transmitter  and  two  reeeivers,  or  multiple  reeeivers  and  multiple 
transmitters.  It  is  a  generalization  of  the  bistatie  radar  system,  with  one  or  more  reeeivers 
proeessing  echo  field  returns  from  one  or  more  geographically  separated  transmitters. 
The  multistatic  radar  system  is  desirable  as  it  allows  covert  operation  of  the  receivers  and 
there  is  increased  resilience  to  electronic  countermeasures.  Additionally,  due  to 
geometrical  effects,  the  radar  cross-section  of  the  target  could  potentially  be  enhanced. 

Conventional  synthetic  aperture  radar  or  inverse  synthetic  aperture  radar  imaging 
involves  the  processing  of  short  pulse  signals  to  form  high  range  resolution  images  of  the 
target.  Similarly,  high  Doppler  resolution  data  could  be  utilized  to  generate  images.  One 
major  benefit  of  employing  high  Doppler  resolution  imaging  in  a  multistatic  radar  system 
is  the  requirement  for  relatively  simple  and  inexpensive  continuous  wave  (CW) 
transmitters  and  receivers.  Eor  a  typical  airborne  target,  the  signal  returns  are 
intrinsically  considered  narrowband  since  the  Doppler  shifts  are  usually  in  the  region  of 
tens  to  hundreds  of  kilohertz.  This  narrowband  attribute  in  a  Doppler  imaging  system  has 
utility  in  operations  where  the  available  frequency  spectrum  is  limited. 

Conceptually,  a  Doppler-only  multistatic  radar  system  could  potentially  provide 
high-resolution  imaging  of  airborne  targets  in  support  of  non-cooperative  target 
recognition  [2],  while  keeping  the  cost  and  ease  of  implementation  and  maintenance  to  a 
significantly  lower  level  when  compared  with  current  pulse-Doppler  monostatic  radars. 
This  is  the  subject  of  this  thesis. 
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B,  OBJECTIVE 

The  objective  of  this  thesis  is  to  develop  a  basic  theory  for  a  Doppler-only 
multistatic  radar  system  based  on  high-Doppler  resolution  data.  The  physics  involved  in 
resolving  the  geometry  of  the  targets  relative  to  the  transmitters  and  receivers,  as  well  as 
the  imaging  process  will  be  discussed.  The  mathematical  implementation  of  such  a  radar 
system  will  also  be  examined.  Thereafter,  the  results  of  the  above  analysis  will  be 
exercised  on  test  cases. 

C.  THESIS  ORGANIZATION 

The  thesis  is  divided  into  eight  chapters,  organized  as  follows. 

Chapter  II  provides  a  review  of  current  imaging  methods,  including  one¬ 
dimensional  high  range  resolution  imaging,  its  extension  to  two-dimensional  imaging, 
and  high  Doppler  resolution  imaging. 

Chapter  III  provides  some  background  of  inverse  problems  in  radar  imaging.  The 
well-posedness  and  condition  number  of  the  problem,  and  the  method  of  least  squares  are 
discussed. 

Chapter  IV  introduces  the  concepts  of  bistatic  radars.  It  presents  an  overview  of 
bistatic  radar  definitions  and  the  parameters  involved  in  the  derivation  of  the  model. 
Thereafter,  it  extends  the  bistatic  model  into  the  multistatic  case.  The  advantages  and 
disadvantages  of  the  multistatic  radar,  as  well  as  the  implementation  requirements,  are 
discussed. 

Chapter  V  introduces  the  design  of  the  radar  model.  It  describes  the  mathematical 
approach  to  derive  the  target’s  position  and  velocity. 

Chapter  VI  involves  the  implementation  of  the  multistatic  radar  model  for 
imaging.  It  illustrates  the  results  and  tabulation  of  measurement  errors. 
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Chapter  VII  provides  the  results  analysis  and  an  examination  of  the  imaging 
artifacts. 

Finally,  Chapter  VIII  recommends  areas  of  future  work  and  concludes  with  some 
comments  on  the  impetus  towards  developing  a  Doppler-only  multistatic  radar  system. 
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II.  RADAR  IMAGING  METHODS 


A.  OVERVIEW  OF  IMAGING  RADAR 

Radar  systems  measure  the  strength  of  the  baekseatter  and  the  round  trip  delay 
time  of  radio  frequeney  signals  refleeted  from  distant  objeets.  Sinee  the  radar  pulse 
travels  at  the  speed  of  light,  it  is  relatively  straightforward  to  use  the  measured  time  delay 
for  the  round  trip  of  a  partieular  pulse  to  calculate  the  range  to  the  reflecting  object.  The 
resolution  in  the  down-range  direction  is  governed  by  the  pulse  bandwidth;  a  higher 
bandwidth  will  imply  finer  resolution  in  this  direction.  The  resolution  of  the  image  in  the 
cross-range  direction  is  determined  by  the  dimensions  of  the  radar  antenna  and  a  larger 
antenna  will  lead  to  finer  resolution  in  the  azimuth  direction. 

This  chapter  discusses  the  current  imaging  methods  such  as  one-dimensional  high 
range  resolution  (HRR)  imaging,  two-dimensional  imaging  and  high  Doppler  resolution 
(HDR)  imaging  will  be  discussed. 

B,  ONE-DIMENSIONAL  IMAGING 

One-dimensional  imaging  involves  the  generation  of  range  profiles,  which  is  a 
straightforward  method  of  articulating  target  substructure.  High  range  resolution  imaging 
involves  the  resolution  of  the  individual  target  scatterers  through  the  down-range  profile 
of  the  target.  This  method  essentially  provides  a  one-dimensional  image  of  the  target  and 
an  example  is  shown  in  Figure  1 .  Short  pulse  or  pulse-compressed  signals  can  be  used  to 
form  one-dimensional  HRR  images  of  the  target.  As  the  short  pulse  sweeps  across  the 
target,  it  sequentially  excites  the  scatterers  on  the  target,  which  re-radiate  energy  back  to 
the  radar  receiver.  When  these  scatterers  are  non-interacting  and  point-like,  the  scattered 
pulse  will  be  a  sum  of  damped  and  blurred  images  of  the  incident  pulse,  which  are  shifted 
by  time  delays  that  are  proportional  to  the  scatterer’s  range  [3].  For  pulse  compression, 
the  two  common  waveforms  are  linear  frequency  modulation  and  binary  phase  coded 
pulses  [4]. 
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Figure  1.  Example  of  a  range  profile  of  a  Boeing  737-500  (From  [5]) 

The  down-range  profile  can  be  affected  by  a  variety  of  factors  such  as  target 
aspect  angle,  position  of  the  scatterers  or  masking  of  scatterers  by  other  parts  of  the 
target.  Additionally,  since  the  bandwidth  of  a  pulse  is  inversely  proportional  to  its 
duration,  the  use  of  short  pulses  will  lead  to  large  bandwidth  requirements.  Wide 
bandwidth  can  increase  system  complexity  and  increase  the  likelihood  of  interference 
from  other  emitters  in  the  electromagnetic  spectrum.  A  short-pulse  waveform  also 
provides  less  accurate  radial  velocity  measurement  as  compared  to  that  obtained  from 
Doppler  frequency  shift.  As  high  peak  power  is  required  to  transmit  short  pulses  over 
long  ranges,  it  is  an  important  limitation  in  radar  applications.  High  peak  power 
transmission  can  result  in  voltage  breakdown,  especially  at  high  frequencies  since  the 
waveguide  dimensions  are  small. 

Radar  imaging  through  using  only  range  profiles  has  limited  applications.  This  is 
due  to  the  fact  a  range  profile  will  not  be  able  to  distinguish  cross-range  target  structures. 
All  scatterers  located  at  the  same  distance  from  the  radar  will  reflect  energy  back  to  the 
radar  with  the  same  time  delay.  Hence,  when  the  radar  illuminates  many  distinct  targets 
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at  any  instant,  conclusive  decisions  pertaining  to  the  nature  of  the  target  cannot  be 
formed  using  a  single  set  of  range-only  data  [6]. 


C.  TWO-DIMENSIONAL  IMAGING 


For  sufficient  target  interpretation,  other  information  is  usually  needed  in  addition 
to  the  range  profile.  Such  added  information  can  be  in  the  form  of  a  high-resolution 
cross-range  profile,  Doppler  profile,  or  the  “triangulation”  of  range  profiles.  The 
Synthetic  Aperture  Radar  (SAR)  and  Inverse  Synthetic  Aperture  Radar  (ISAR)  are 
examples  of  two  different  sehemes  of  imaging  systems  that  form  two-dimensional  images 
of  a  target. 

“Triangulation”  of  different  sets  of  range  profiles  allows  the  extension  of  the  radar 
imaging  concept  to  two  dimensions.  This  approach  permits  the  determination  of  cross¬ 
range  target  structure  while  using  only  HRR  radar  systems.  It  relies  on  collecting 
multiple  sets  of  range  profiles  from  different  target  orientations,  processing  them,  and 
synthesizing  an  image. 

Figures  2  and  3  illustrate  this  eoncept.  Consider  a  set  of  three  point  targets  with 
the  radar  located  at  the  same  distance  from  targets  2  and  3  (Figure  2).  When  the 
geometry  is  as  in  Figure  2,  the  return  echo  will  only  indicate  two  targets.  This  illustrates 
the  eondition  where  ambiguity  exists  when  targets  lie  along  the  bands  of  constant  range 
from  the  radar. 


Target  1 


Target  2 


Target  3 


Bands  of  constant 
range  from  radar 


Figure  2. 


Echo  signal  shows  only 
two  targets  in  range 


Ambiguity  scenario  from  a  single  pulse 
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With  multiple  sets  of  data  from  different  direetions,  triangulation  will  allow  for 
gradual  buildup  of  the  relative  positions  of  the  three  targets  (Figure  3).  The  range 
profiles  are  swept  in  the  cross-range  direction  to  form  bands  of  constant  range  from  the 
radar,  which  represent  the  possible  locations  of  the  target  scatterer.  These  bands  are  then 
superimposed  and  the  crossing  points  are  used  to  determine  the  scattering  center 
locations.  With  correct  correlation  methods,  a  target  image  built  up  of  points  where  the 
swept  lines  intersect  can  be  obtained. 


Figure  3.  Cross-range  information  obtained  from  range  profiles 

In  essence,  two  different  schemes  are  used  to  collect  target  data  from  different 
target  aspects  to  obtain  cross-range  information.  The  radar  can  move  around  the  target 
while  it  remains  fixed;  or  the  target  can  rotate  while  the  radar  stares  at  it  over  time.  In  the 
former  case,  the  radar  is  described  as  collecting  data  over  a  synthetic  aperture,  and  thus, 
termed  Synthetic  Aperture  Radar  (SAR).  In  the  latter  case,  the  radar  is  described  as  an 
Inverse  Synthetic  Aperture  Radar  (ISAR).  Mathematically,  the  basics  for  image 
formation  are  fundamentally  the  same  for  both  schemes.  Notwithstanding,  the  details  of 
data  processing,  such  as  correcting  for  departures  from  ideal  target  behavior,  are  different 
in  the  two  settings. 
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D.  HIGH  DOPPLER  RESOLUTION  (HDR)  IMAGING 


The  concept  of  utilizing  high  Doppler  resolution  inputs  to  synthesize  a  target 
image  can  be  illustrated  by  drawing  an  analogy  from  HRR  imaging. 

In  HRR  imaging,  the  image  is  developed  by  a  correlation  of  range  profiles  taken 
at  different  locations.  Each  “line”  of  the  range  profile  is  an  iso-range  line  at  a  specific 
time  since  the  radar  pulse  takes  a  finite  amount  of  time  to  reflect  off  a  scatterer  and  return 
to  the  receiver. 

In  HDR,  another  set  of  “lines,”  the  iso-Doppler,  which  are  lines  of  constant 
Doppler  shift  at  a  specific  frequency,  is  used  to  build  the  image.  For  a  target  rotating 
about  its  center,  the  local  velocity  of  a  scatterer  will  depend  on  its  radial  distance  from 
the  center.  The  return  signal  at  a  given  Doppler  shift  is  a  superposition  of  all  the  returns 
due  to  scatterers  with  the  same  closing  velocity  lying  along  the  iso-Doppler  hyperbola. 
Conceptually,  the  iso-range  lines  and  iso-Doppler  lines  can  be  considered  as  taking 
“slices”  of  the  image,  with  the  difference  between  the  two  being  the  “slices”  are 
orthogonal  to  each  other,  as  illustrated  in  Figure  4  [2]. 


HRR  down-range  “slices” 
at  time  intervals 


HDR  cross-range  “slices” 
at  frequency  intervals 


Figure  4.  Orthogonal  “slices”  of  the  target  between  using  HRR  and  HDR  (After  [2]) 
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Similarly,  ambiguities  due  to  multiple  seatterers  along  an  iso-Doppler  ean  be 
resolved  by  aequiring  multiple  Doppler  profdes  from  different  aspeet  angles.  However,  it 
should  be  noted  that  the  HDR  eross-range  sliees  only  works  for  rotating  targets. 

In  HRR  imaging,  eaeh  range  profile  has  high  resolution  in  the  down-range 
direetion  and  the  eorrelation  of  the  profiles  for  different  aspeet  angles  will  provide  the 
eross-range  resolution  to  build  a  two-dimensional  image.  Conversely,  in  HDR  imaging, 
eaeh  Doppler  profile  will  have  high  resolution  in  the  eross-range  direetion  and  the 
eorrelation  proeess  will  provide  the  down-range  resolution. 
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III.  INVERSE  PROBLEMS  IN  RADAR  IMAGING 


A.  GENERAL  FORMULATION  OF  THE  IMAGING  PROBLEM 

A  general  formulation  of  the  imaging  problem  ean  be  achieved  by  considering  a 
space-invariant  imaging  system  represented  by  the  operator  S,  with  input  image  f(x)  and 
output  image  g(x).  They  are  related  by  Equation  (3.1)  and  the  model  is  illustrated 
schematically  in  Figure  5.  Since  the  system  S  is  known,  a  real  non-negative  image  g(x) 
with  values  in  a  given  domain  x  e  91 ,  is  generated  at  the  output  of  the  system. 

g  =  Sf  (3.1) 


> 
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Figure  5.  Image  formation  model 


The  image  synthesis  problem  first  determines  if  an  image  /  e  3  (that  satisfies  the 
given  constraints  represented  by  the  set  3)  exist  such  that  the  system  output  Sf(x)  =  g(x) 
for  all  X  e91  [7].  If  such  an  image  does  exist,  then  the  next  step  is  to  determine  f.  On  the 
other  hand,  if  the  image  does  not  exist,  then  the  approach  is  to  determine  an  approximate 
input  image  /  that  produces  an  output  g  =  Sf  according  to  some  criterion.  One  possible 
optimization  criterion  for  finding  an  approximate  solution  is  to  minimize  a  norm 

The  solution  to  this  general  formulation  will  “define”  the  solution  of  inverse  problems. 

B,  INVERSE  PROBLEMS 

J.B.  Keller  formulated  the  following  general  definition  [8]  of  inverse  problems, 
which  is  often  cited  in  literature: 
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We  call  two  problems  inverses  of  one  another  if  the  formulation  of  each 
involves  all  or  part  of  the  solution  of  the  other.  Often,  for  historical 
reasons,  one  of  the  two  problems  has  been  studied  extensively  for  some 
time,  while  the  other  is  newer  and  not  so  well  understood.  In  such  cases, 
the  former  problem  is  called  the  direct  problem,  while  the  latter  is  called 
the  inverse  problem. 

A  direct  problem  can  be  thought  of  as  one  oriented  towards  a  cause  and  effect 
sequence,  while  the  corresponding  inverse  problem  is  linked  to  the  reversal  of  this  cause 
and  effect  sequence  and  comprises  determining  the  unknown  causes  of  known 
consequences.  In  radar  imaging,  the  direct  problem  is  the  mapping  from  the  target  to  the 
quantities  that  can  be  measured  by  the  radar. 

Consequently,  the  inverse  problem  is  the  problem  of  finding  the  original  target 
given  the  data  and  knowledge  of  the  direct  problem.  This  mapping  is  called  the  target 
image.  However,  this  problem  is  ill-posed  due  to  the  loss  of  information  intrinsic  to  the 
solution  of  the  direct  problem.  Therefore,  if  an  image  corresponds  to  two  distinct  objects, 
then  the  solution  of  the  inverse  problem  is  not  unique.  Further,  if  there  are  two 
neighboring  images  such  that  the  corresponding  objects  are  very  distant,  then  the  solution 
of  the  inverse  problem  does  not  depend  continuously  on  the  data. 

The  accepted  approach  for  solving  inverse  problems  which  are  ill-posed  is  to 
search  for  approximate  solutions  satisfying  additional  constraints  based  on  the  physics  of 
the  problem  [12].  This  set  of  approximate  solutions  corresponding  to  the  same  data 
function  is  the  set  of  objects  with  images  close  to  the  measured  one. 

C.  WELL-POSED  AND  ILL-POSED  PROBLEMS 

The  mathematical  term  for  well-posed  problems  arose  from  a  classical  concept 
defined  by  Hadamard.  He  believed  that  such  problems  are  considered  well-posed  if  they 
satisfy  the  following  three  conditions  [10]; 

1 .  the  solution  is  unique 

2.  the  solution  exists  for  any  data 

3 .  the  solution  depends  continuously  on  the  data  and  parameters 
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If  one  of  the  three  conditions  above  is  not  satisfied,  the  problem  is  said  to  be  ill- 
posed.  Hence,  an  ill-posed  problem  is  one  whose  solution  is  not  unique  and/or  does  not 
exist  for  any  data  and/or  does  not  depend  continuously  on  the  data. 

The  third  condition  is  motivated  by  the  fact  that  the  mere  process  of  measuring 
data  involves  small  errors.  If  the  problem  is  well-posed,  the  error  propagation  is 
controlled  by  the  condition  number.  Consider  the  solution  of  a  linear  system  of  algebraic 
equations 

Ax  =  m  (3.2) 

where  m  is  the  measured  data,  x  is  the  signal  and  A  is  an  operator  that  describes  the 
nature  of  the  measurement  system.  In  the  case  of  a  unique  solution,  if  there  are  n 
variables,  then  there  must  be  n  independent  equations  in  order  to  obtain  a  solution.  By 
conditions  1  and  2,  the  inverse  of^  exists.  By  thinking  of  Am  as  being  the  small  error  in 
the  measured  data  and  Ax  as  being  the  resulting  error  in  x  such  that 

A(x  +  Ax)  =  m  +  Am 

Consequently,  this  uncertainty  Ax  complies  with 

||Ax||  ||A»i|| 

-j— j-  <  cond{A)— — ^ 

IfII  II"*II 

where  ||Am||/||m||  is  the  relative  change  in  measurement  and  ||Ax||/||x||  is  the  relative 
error  caused  by  this  change  [11]. 

The  condition  number  is  a  relative  error  magnification  factor  and  the  quantity  that 
controls  error  propagation  from  data  to  solution.  It  is  defined  as 

cond(A)  =  11^1 1^^' I 

where  the  common  norm  of  Euclidean  distance  is  adopted  for  ||^||  such  that 


1 

2 
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It  should  be  noted  that  the  computation  of  the  matrix  norm  corresponding  to  the  /2-norm 
criterion  involves  the  singular  value  decomposition  (SVD).  If  the  condition  number  of  A 
is  not  large,  then  the  problem  is  considered  to  be  well-conditioned.  The  solution  to  the 
problem  will  be  stable  relative  to  small  variations  of  data.  This  requirement  of  stability  is 
essential  for  meaningful  problems  in  approximation  methods.  On  the  other  hand,  if  the 
condition  number  is  large,  then  the  problem  is  deemed  ill-conditioned.  The  ill- 
conditioned  situation  is  somewhat  similar  to  the  ill-posed  condition.  An  ill-posed 
problem  can  be  thought  of  as  one  where  an  inverse  does  not  exist  because  m  +  Am  is 
outside  the  range  of  A. 

In  general,  most  inverse  problems  are  ill-posed.  In  the  case  of  a  bandlimited 
system,  the  solution  of  the  inverse  problem  is  not  unique  and  the  first  condition  required 
for  well-posedness  is  not  satisfied.  This  is  because  the  imaging  system  does  not  transmit 
information  about  the  target  at  frequencies  outside  the  band  of  the  instrument. 

Straightforward  solutions  of  ill-posed  problems  can  result  in  non-physical 
answers.  It  is  interesting  to  note  that  a  more  refined  operator  A  may  work  against 
obtaining  a  more  reliable  solution.  When  discretizing  an  ill-posed  problem,  the  condition 
number  of  the  corresponding  discrete  problem  can  be  very  large.  It  is  the  case  that,  when 
the  discretization  of  the  ill-posed  problem  gets  finer,  the  condition  number  of  the 
corresponding  discrete  problem  becomes  larger  [9].  For  this  reason,  the  discretization  of 
equations  must  be  done  carefully  and  techniques  must  be  introduced  to  construct  stable 
solutions. 

D,  METHOD  OF  LEAST  SQUARES 

An  overdetermined  system  is  one  that  comprises  more  equations  than  unknowns. 
Overdetermined  systems  of  simultaneous  linear  equations  are  often  encountered  in 
various  kinds  of  curve  fitting  to  experimental  data.  The  method  of  least  squares  is 
frequently  used  to  solve  such  systems  of  equations  in  an  approximate  sense.  During  the 
process  of  measuring  data,  errors  or  inaccuracies  will  be  manifested  in  the  measurement. 
Hence,  instead  of  solving  the  equations  accurately,  the  minimization  of  the  sum  of 
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squares  of  the  residuals  is  sought.  The  result  is  a  eurve  fit  to  a  set  of  data  points  sueh  that 
the  sum  of  squares  of  the  range  between  the  data  point  and  its  eorresponding  eurve  is 
minimized. 

If  appropriate  probabilistie  assumptions  about  underlying  error  distributions  are 
made,  the  least  squares  approaeh  beeomes  what  is  known  as  the  maximum-likelihood 
estimate  of  the  parameters.  This  least  squares  eriterion  (f  -norm  eriterion)  is  widely  used 
for  the  resolution  of  inverse  problems  due  to  its  simplieity.  It  should  be  noted  that  at 
times,  it  is  applied  even  if  its  basie  underlying  hypothesis  (Gaussian  uneertainties)  is  not 
always  satisfied  [13]. 

In  the  eurve  fitting  model,  let  t  be  the  independent  variable  and  let  m(t)  denote  the 
unknown  funetion  of  t  to  be  approximated.  If  there  are  k  observations,  then  the  values  of 
m  measured  at  speeified  points  of  t  are 

mi  =  m(ti),  i  =  \,2,...,k 

The  funetion  m(t)  is  then  modeled  by  a  linear  eombination  of  n  basis  funetions 
m(t)  Xi<p^(t) +  ....  +  xj„(t) 

The  design  matrix  A  is  a  reetangular  matrix  of  order  khy  n  with  elements 

ay  = 

In  matrix-veetor  notation,  the  model  is 
m  «  Ax 

There  are  many  models  available  for  the  eurve  fitting  solution.  Some  eommon  models 
inelude  the  straight  line,  polynomials,  rational  funetions,  exponentials,  log-linear  and 
Gaussian.  This  thesis  will  foeus  on  applying  the  polynomial  model  for  the  least  squares 
applieation.  The  model  used  is  a  eubie  polynomial  of  the  form 

m{t)  «  Xj  +  Xjt  +  xd^  +  xd^ 

The  residuals  are  the  differenees  between  the  measurements  and  the  model; 
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j=i 

=  m.  -  m{t) 

=  m.  -  +  x^t  +  x^t^  +  x^t^ ^ 

By  the  method  of  least  squares,  the  objeetive  is  to  minimize  the  sum  of  squares  of  the 
residuals: 

k 

II  if  V  2 

IHI  =Z.f 

i=\ 

k  2 

^  Ikf  “  “  (-^1  +-^2^  +  -^3^^  +-^4^^)]  (3-3) 

/=! 

E,  CORRELATION  RECEPTION 

In  radar  processing,  correlation  reception  is  an  important  technique  that  allows  the 
accurate  separation  of  the  desired  signal  from  the  unwanted  noise.  It  involves  the 
comparison  of  the  received  signal  with  reference  signals  of  the  form 

=  (3.4) 

where  Sscatt(t’)  is  the  radar  signal  associated  with  the  field  reflected  from  an  object  at 
range  t  and  radial  velocity  v  and  /?  is  a  signal  strength  scale  factor  [14],  This  comparison 
is  performed  coherently  such  that  the  phase  of  the  transmission  signal  is  preserved  in  the 
reference  signal.  The  correlation  receiver  can  be  considered  as  a  type  of  matched  filter 
receiver. 

The  original  radar  signal  processing  problem  involves  optimal  detection  in 
additive  noise  so  that  the  received  signal  is  of  the  form 

Sreci^)  =  S,,,„{^)+n{t)  (3.5) 

where  n(t)  is  a  random  noise  process  and  Srec(t)  is  the  received  signal  measured  by  the 
radar.  The  objective  is  to  determine  Sscatt(t)  from  the  random  measurements  of  Srec(t)- 
Estimation  is  usually  accomplished  through  maximum  likelihood  processing  whereby 

Srec(t)  is  compared  to  some  idealized  signals  generated  using  a  signal  model.  The 
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maximum  likelihood  estimate  for  Sscatt(t)  is  one  which  maximizes  the  probability  , 
which  is  the  conditional  probability  density  function  of  the  measured  signal  for  a  specihc 

scattered  field  [6]. 


e 

^  scan,  ML 


=  arg  max  p  ,  (u) 

u  G  model  space  ^  I 


Consequently,  it  can  be  shown  that  correlation  receivers  attempt  to  determine  the  time 
shift,  T ,  and  frequency  shift,  o  ,  which  maximizes  the  real  part  of  the  correlation  integral 

r/{u,T)=\  (3.6) 

J  —00 


In  radar  operations,  the  natural  model  is  based  on  the  scattering  interaction 
between  the  interrogating  held  and  target.  If  Si„c(t)  is  the  incident  pulse  transmitted  by 
the  radar,  then  the  linear  radar  scattering  model  is  derived  from  the  superposition  of 
Equation  (3.4). 

^scaAt)  =  f  f”  -  T)e“^^‘^^'>dTdu  (3.7) 

—00 


Now,  this  model  is  substituted  into  Equation  (3.5)  and  then  into  the  correlation  integral  of 
Equation  (3.6). 


r/io,  r)  =  j  j  T)slXt'-  T')e 


dt'dr'du' 


+  correlation  noise  term  (3.8) 

The  correlation  noise  term  is  considered  small  and  can  thus  be  ignored.  Thereafter, 
Equation  (3.8)  can  be  rearranged  to 


p{v\T')x{v-v',T-T')e 

J  J  —00 


i  —  {v+v'){T-T') 

2  dr'dv' 


(3.9) 


where 


X(u,t)  =  ^  s.. 


(3.10) 
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Equation  (3.9)  is  the  standard  radar  data  model  and  expresses  the  output  of  the  correlation 
receiver  as  the  convolution  of  p  and  x  ■  The  function  x(^,  defined  by  Equation  (3.10) 
is  the  ambiguity  function. 

F.  AMBIGUITY  FUNCTION 

The  accuracy  with  which  the  target  position  and  radial  velocity  can  be  estimated 
from  the  radar  data  is  determined  by  the  ambiguity  function.  In  radar  systems, 
waveforms  are  selected  to  optimize  the  requirements  of  detection,  measurement  accuracy, 
resolution,  and  ambiguity  and  clutter  rejection.  Ambiguity  function  plots  are  scrutinized 
for  a  qualitative  determination  of  the  suitability  of  different  waveforms  in  meeting  the 
above  requirements  [4].  In  practice,  the  ambiguity  function  is  plotted  as  a  function  of 
time  delay  and  Doppler  shift.  In  radar  imaging,  the  ambiguity  function  can  be  thought  of 
as  the  imaging  kernel. 

1.  Ambiguity  Function  for  HRR  Signal 

In  HRR  imaging,  the  system  can  be  represented  by  an  impulse  radar  transmitting 
a  very  short  pulse  in  the  form  of  a  Dirac  delta  function 

=  S{t) 

The  ambiguity  function  becomes 

X{v,  r)  =  J  -  r)  5,„^(t'+  -  T)e‘^‘  dt' 

2  2 

/•  00  T  T  , 

=  f  df 

J-*  2  2 

.  r 

.-.  x{<^,  T)  =  d(T)e“'^ 

The  data  collected  from  such  a  radar  is  independent  of  frequency  shift  and  is 
known  as  the  down-range  profile.  In  reality,  a  Dirac  delta  function  does  not  exist 
physically,  hence,  the  waveform  used  is  usually  an  approximation  of  an  impulse. 
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2,  Ambiguity  Function  for  HDR  Signal 

In  HDR  imaging,  the  system  ean  be  represented  by  a  eontinuous  wave  radar 
transmitting  a  single  tone  signal  in  the  form 

The  ambiguity  funetion  then  beeomes 

•'-*  2  2 
=  \  e  ^  e"^‘ dt' 

J  —00 

=  e^‘dt' 

J  —00 

=  (3.11) 

It  ean  be  eoneluded  that  the  data  eolleeted  is  independent  of  down-range  delay. 
Consequently,  sueh  radar  data  ean  be  used  to  determine  the  distribution  of  down-range 
veloeity  aeeurately,  but  will  offer  no  information  about  target  position. 

3.  Bistatic  Ambiguity  Function 

For  a  monostatic  radar,  the  use  of  delay  and  Doppler  shift  as  the  arguments  of  the 
ambiguity  function  is  sufficient  due  to  the  linear  relationship  between  target  range  and 
range  rate.  Thus,  the  delay  and  Doppler  shift  pair,  or  the  range  and  range  rate  pair,  may 
be  used  interchangeably  in  the  ambiguity  function. 

In  the  case  of  a  bistatic  or  multistatic  radar,  the  situation  is  more  complicated. 
Since  the  transmitter  and  receiver  are  not  at  the  same  site,  the  relationship  between 
Doppler  shift  and  target  velocity,  and  between  time  delay  and  range,  are  highly  non¬ 
linear.  The  shape  of  the  ambiguity  function  is  a  strong  function  of  the  geometry  between 
the  transmitter,  receiver  and  target,  as  well  as  waveform  properties  [15].  Hence,  the 
representation  of  the  ambiguity  function  for  a  bistatic  radar  in  terms  of  delay  and  Doppler 
shift  is  not  meaningful  and  may  even  lead  to  incorrect  conclusions. 
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It  has  been  proposed  that  the  ambiguity  funetion  for  bistatie  radar  be  written  as: 


XiRRH^RRa^V^j,V^,e^,L)  = 

00 

j  s{t-TXRRa,OR,L)-s\t 


1 


[-  i  {cojyH  ( ^RH  fi, 6 n.L)- cojy^  {Rji^ ,  ,  Z, ))?  ] 


dt 


(3.12) 


where  Rrh  and  RRa  are  the  hypothesized  and  aetual  ranges  from  the  reeeiver  to  the  target; 
Vh  and  Va  are  the  hypothesized  and  aetual  radial  veloeities  of  the  target  relative  to  the 
receiver;  and  are  the  hypothesized  and  actual  Doppler  shift;  6^  is  the  bearing 

of  the  target  with  respect  to  the  receiver;  and  L  is  the  baseline  length  between  the 
transmitter  and  receiver  [16].  Equation  (3.12)  assumes  the  reference  point  of  the  bistatic 
geometry  to  be  the  receiver.  The  main  difference  between  Equations  (3.12)  and  (3.10)  is 
that  the  geometrical  disposition  of  the  transmitter,  receiver,  and  target  are  taken  into 
account. 
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IV.  BISTATIC  AND  MULTSTATIC  RADAR  THEORY 


A.  BACKGROUND 

In  monostatic  radar  systems,  increased  information  is  frequently  associated  with 
increased  bandwidth.  With  the  increased  bandwidth,  high  range  resolution  is  obtainable 
and  frequency  diversity  offers  one  such  approach  to  obtain  additional  information  about 
targets.  On  the  other  hand,  in  a  multistatic  radar  system,  geometric  diversity  [17]  in  the 
disposition  of  transmitting  and  receiving  stations  allows  additional  target  information  to 
be  obtained.  With  suitable  signal  processing,  the  additional  information  obtained  can 
translate  into  improved  target  detection.  Additionally,  geometric  diversity  offers  the 
potential  for  increased  resolution  and  it  is  a  dual  to  frequency  diversity  in  classical 
monostatic  radar  systems.  Hence,  improved  target  position  accuracy  and  image 
resolution  can  be  expected  in  a  multistatic  radar  system. 

The  image  resolution  for  the  monostatic  synthetic  aperture  radar  (SAR)  can  be 
drawn  from  tomography/Fourier  space  or  radar  range/Doppler  principles.  In  the  case  of 
multistatic  SAR,  the  resolution  result  is  understood  using  the  tomography/Fourier  space 
principles.  As  an  example,  an  ultra  wideband  monostatic  SAR  with  a  50%  bandwidth 
will  have  a  range  resolution  of  X  and  a  pixel  area  of  .  A  Multistatic  SAR  system  will 
have  a  range  resolution  of  X  !  3  and  a  pixel  area  of  X^  !  9  [18].  This  translates  to  a  9.5 
dB  improvement  over  monostatic.  A  similar  performance  gain  can  be  expected  from  a 
multistatic  ISAR  system. 

Since  the  bistatic  radar  is  the  building  block  for  the  multistatic  radar  system,  the 
fundamental  theory  for  bistatic  radar  operation  will  first  be  discussed.  These  will  include 
the  coordinate  system,  the  determination  of  target  location  and  Doppler  information  and 
resolution.  Thereafter,  the  extension  of  the  bistatic  radar  to  the  multistatic  radar  system, 
its  advantages  and  drawbacks,  and  implementation  requirements  will  be  examined. 
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B. 


DEFINITION  OF  BISTATIC  RADAR 


As  defined  by  the  IEEE  Standard  686-1997  [1],  the  bistatic  radar  is  a  radar  using 
antennas  for  transmission  and  reception  at  sufficiently  different  locations  that  the  angles 
or  ranges  from  those  locations  to  the  target  are  significantly  different.  Bistatic  radar  is 
deemed  to  have  significant  advantages  against  stealthy  targets  due  to  the  relatively  larger 
amount  of  forward  scattering.  Furthermore,  the  operation  of  the  bistatic  radar  receiver  is 
considered  covert,  as  its  location  cannot  be  determined  easily  by  merely  observing  the 
radar’s  emissions. 

Figure  6  illustrates  the  setup  of  a  bistatic  radar,  with  a  transmitter  Tx  and  receiver 
Rx  situated  at  two  different  locations,  separated  by  a  baseline  distance  L.  The  bistatic 
angle  P  is  the  angle  between  the  transmitter  and  receiver  with  the  vertex  at  the  target. 


Target,  Tgt 


Transmitter,  Tx 


Reeeiver,  Rx 


Figure  6.  Bistatic  radar 

C.  COORDINATE  SYSTEM 


This  thesis  will  adopt  a  two-dimensional.  North-referenced  coordinate  system  for 
derivation  and  discussion.  Figure  7  depicts  the  coordinate  system  and  the  parameters 
characterizing  bistatic  radar  operation.  The  plane  containing  the  transmitter,  receiver  and 
target  is  defined  as  the  bistatic  plane  and  the  bistatic  triangle  lies  in  this  plane.  The 
angles,  <l)j  and  are  the  respective  target  bearings  for  the  transmitter  and  receiver. 
They  are  taken  to  be  positive  when  measured  clockwise  from  North.  The  angle  /  2  is 
known  as  the  bistatic  bisector  angle.  It  should  be  noted  that  . 
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Figure  7.  Bistatic  radar  (North  referenced)  coordinate  system 


D.  TARGET  LOCATION 

The  target’s  position  relative  to  the  receiving  station  can  be  determined  if  Rr  and 
can  be  obtained.  However,  the  receiver-to-target  range,  Rr,  cannot  be  measured 
directly  as  a  bistatic  radar  usually  measures  target  range  as  the  range  sum  Rt  +  Rr.  It  can 
be  calculated  by  solving  the  bistatic  triangle  of  Figure  7.  From  the  cosine  rule,  Rt  and  Rr 
are  calculated  as  follows: 


Rj  —  Rjf  -\-  L  —  2,RjfL  cos|^~  + 

RI  +  2Rj^L  sin ((pj^ )=  R^  - 

2Rj^  +  2Rj^L  sin^^^^+  27?^  7?^  =  (Rj.  +  —  L 


,  (r,  +  rJ-l^ 

2[(7?^+7?^)+Lsin(^^)] 


(4.1) 


From  equation  (4.1),  Rr  can  be  determined  once  the  range  sum  Rr  +  Rr  is 
obtained  and  (pj^  is  known.  Thereafter,  Rr  can  be  expressed  as  follows: 


Rj^  =  ^JrI  +L^  -  2Rj^L  sin(^zi^  ) 


23 


(4.2) 


Alternatively,  Rj  can  also  be  expressed  in  terms  of  the  range  sum  Rj  +  Rr  and  (j)j 


by  making  use  of  the  eosine  rule. 

^  {Rt+RrI-l^ 


(4.3) 


The  range  sum  Rr  +  Rr  ean  be  determined  by  two  methods,  namely  direet  and 
indireet.  In  the  direct  method,  the  reeeiver  measures  the  time  interval  AT),,  between 
reeeption  of  the  transmitted  pulse  and  reeeption  of  the  target  eeho.  Adequate  line  of  sight 
between  transmitter  and  reeeiver  is  required.  The  range  sum  is  then  ealeulated  as 

{Rt+R^)=cM^  +  L  (4.4) 

In  the  indirect  method,  synehronized  stable  eloeks  are  used  by  the  reeeiver  and  a 
dedieated  transmitter.  Line  of  sight  between  transmitter  and  reeeiver  is  not  required  for 
this  measurement.  The  reeeiver  measures  the  time  interval  AT)  between  transmission  of 
the  pulse  and  reeeption  of  the  target  eeho.  The  range  sum  is  then  ealeulated  as 

{RT+R^)=c^TJ  (4.5) 


For  this  thesis,  the  direct  method  of  time  interval  measurement  is  adopted. 

Additionally,  from  the  sine  rule,  the  following  relationships  are  obtained: 


sin(y0) 

L 


^  j. 

sini  -  + 
R^ 


Rj^  _  cos{(pj.) 

Rj  cos{(/)g) 

Rj^  _  cos{(pj.) 
L  sin(y0) 


Rj  _  eos(^zJ^) 
L  sin(y0) 


(4.6) 

(4.7) 

(4.8) 

(4.9) 
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E. 


TARGET  DOPPLER 


The  bistatic  Doppler  shift, is  defined  as  the  time  rate  of  change  of  the  total  path 
length  of  the  scattered  signal,  normalized  by  the  wavelength  X .  Since  the  total  path 
length  is  the  range  sum  of  Rj  +  Rr,  the  target  Doppler  shift  is  as  follows. 


X\dt 


dRr  dR, 

■  +  ■ 


X\  dt 


dt 


(4.10) 


Tgt 


Figure  8.  Bistatic  Doppler  geometry 


In  the  case  when  the  transmitter  and  receiver  are  stationary  while  the  target  is 
moving  (see  Figure  8),  the  target’s  bistatic  Doppler  shift  at  the  receiver  station  can  be 
expressed  as  a  function  of  the  aspect  angle  referenced  to  the  bistatic  bisector,  5 .  The 
corresponding  Doppler  shift  is  derived  as  follows; 


dRj. 

dt 


=  Fcos 


(4.11) 


dRj^ 

dt 


=  V  cos 


(4.12) 
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(4.13) 


f=L 


\  ^  { s:  P 

COS  O - +COS  — 


2V 


fB=—cos(S)cos  - 


A 


\2 


The  properties  of  Equation  (4.13)  are  presented  in  Table  1. 


Condition 

Property 

All  p 

Magnitude  of  the  bistatie  target  Doppler  shift  is  never  greater  than 
that  of  the  monostatic  target  Doppler  shift,  when  the  monostatic 
radar  is  located  on  the  bistatie  bisector. 

All  p 

-90°  <  ^  <  +90° 

Bistatie  target  Doppler  shift  is  positive. 

All  p 

5  =  ±90° 

Bistatie  Doppler  shift  is  zero. 

All  <  180° 

^  =0°or  180° 

Target’s  velocity  is  collinear  with  bistatie  bisector.  Magnitude  of 
bistatie  target  Doppler  shift  is  maximum. 

II 

oo 

o 

0 

Target  is  on  the  baseline,  characterized  by  large  forward  scatter. 
Radar  ean  only  deteet  presence  of  a  target,  the  range  is 
indeterminate  due  to  eelipsing  of  the  direet  pulse  by  the  scattered 
pulse. 

Table  1.  Properties  of  bistatie  Doppler  shift  under  various  eonditions  (From  [22]) 


F.  DOPPLER  RESOLUTION 

For  both  monostatie  and  bistatie  radars,  the  Doppler  resolution  depends  on  the 
amount  of  Doppler  separation  between  two  target  echoes  at  the  reeeiver,^^;  and  fsi-  It  is 
accepted  to  be  1/T,  where  T  is  the  receiver’s  coherent  processing  interval  [20].  Henee, 
the  requirement  for  Doppler  resolution  is 

h,-f.2=J  (4-14) 

In  the  bistatie  case,  an  example  of  the  geometry  for  two  targets  is  as  shown  in 
Figure  9.  These  targets  are  assumed  to  be  collocated;  consequently,  they  share  a  common 
bistatie  bisector.  The  corresponding  bistatie  target  Doppler  shifts  are 
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far- 

fB2  = 


2V^ 

2V, 

A 


cos((5'j)cos[^ 

cos(<5’2)cos 


Combining  these  equations,  the  bistatic  Doppler  resolution  can  be  obtained  from 
the  difference  between  the  two  target  velocity  vectors,  projected  onto  the  bistatic  bisector. 

AF  =  fj  cos(Aj)  -  V,  cos((52) 


^AV= 


A 


2Tcos 


(4.15) 


Tgt  #1,  #2 


Figure  9.  Bistatic  Doppler  resolution  (for  2  collocated  targets) 

G.  DEFINITION  OF  MULTISTATIC  RADAR 

The  IEEE  Standard  686-1997  [1]  has  defined  the  multistatic  radar  as  a  radar 
system  having  two  or  more  transmitting  or  receiving  antennas  with  all  antennas  separated 
by  large  distances  when  compared  to  the  antenna  sizes.  It  comprises  at  least  three 
components  -  one  transmitter  and  two  receivers,  or  multiple  receivers  and  multiple 
transmitters.  It  is  an  extension  of  the  bistatic  radar  system,  with  one  or  more  receivers 
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processing  returns  from  one  or  more  geographieally  separated  transmitters.  Examples  of 
sueh  radar  systems  inelude  the  Air  Foree  Spaee  Surveillanee  System  and  the  Silent 
Sentry  System.  Figure  10  depiets  an  illustration  of  a  multistatie  radar  set  up. 


Tgt 


Figure  10.  Sample  Multistatie  radar  eonfiguration  (1  transmitter  &  2  reeeivers) 

H,  MAIN  ADVANTAGES 

Due  to  information  fusion  from  several  spatially  separated  reeeiving  stations,  the 
multistatie  radar  has  signifieant  advantages  over  a  basie  monostatie  radar.  Some  of  these 
advantages  inelude  aeeurate  target  position  estimates,  Doppler  estimation,  eoordinate 
measurement,  inereased  resolution,  power  advantages,  flexible  eonfiguration,  and 
improved  elutter  rejeetion  [21]. 
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1,  Accurate  Position  Estimate  of  Targets 

In  a  conventional  monostatic  radar,  target  position  determination  is  more  accurate 
in  down-range  than  in  eross-range  direetion.  The  multistatie  radar  allows  the  estimation 
of  the  target  eoordinates  through  range-sum  measurements  relative  to  the  spatially 
separated  transmitting  and  reeeiving  stations. 

The  track  update  rate  may  usually  be  higher  in  a  multistatic  radar  than  in  a 
monostatic  radar,  and  consequently,  higher  traeking  accuracy  can  be  aehieved. 

2,  Doppler  Estimation  of  Target  Velocity  and  Acceleration 

Doppler  frequeney  shifts  measured  at  the  various  receiving  stations  allow  the 
estimation  of  the  target’s  velocity  vector.  With  Doppler  information,  trajeetory 
information  of  a  target  ean  be  estimated  with  aceuraey  in  a  short  time  interval.  The  use 
of  Doppler  frequeney  shifts  for  veloeity  and  aeeeleration  estimates  improve  traeking 
aeeuraey  and  general  quality  of  the  traeking  proeess. 

3,  Measurement  of  Three  Dimensional  Coordinates  and  Velocity 

Monostatic  or  bistatic  radars  can  only  determine  the  signals’  Directions  Of 
Arrival  (DOA)  based  on  the  bearings  of  the  radiation  sources.  Multistatic  radars  can 
obtain  the  three-dimensional  coordinates  and  their  derivatives  by  triangulation, 
hyperbolic  methods,  or  their  combination.  Triangulation  determines  the  position  based 
on  the  intersection  of  DOA  from  various  receiving  stations.  The  hyperbolic  method 
determines  the  position  based  on  the  intersection  of  hyperboloids  of  revolution,  which 
have  their  foci  at  receiving  stations.  A  fixed  Time  Difference  Of  Arrival  (TDOA)  at  a 
pair  of  stations  determines  a  hyperboloid  of  revolution  on  which  surface  the  target  must 
lie.  This  TDOA  is  estimated  by  the  signal  delay  in  one  station,  which  is  necessary  to 
maximize  the  mutual  correlation  of  signals  received  by  the  two  stations  [21]. 

The  measurement  of  Doppler  frequency  shifts  of  the  mutual  correlation  function 
of  signals  received  by  a  pair  of  receiving  stations  from  a  moving  target  allows  the 
multistatic  radar  to  estimate  the  radial  velocity  relative  to  these  stations.  A  multistatic 
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radar  system  comprising  four  or  more  stations  will  be  able  to  obtain  all  three  components 
of  the  target  velocity  vector  by  Doppler  frequency  shift  measurements. 

4,  Increased  “Signal  Information”  of  Target  Body 

“Signal  information”  [21]  refers  to  the  information  extracted  from  target  echoes 
that  concerns  geometrical,  physical  and  other  features  such  as  movement  relative  to  target 
center  of  mass.  When  a  target  is  observed  by  a  multistatic  radar  system  from  different 
directions  simultaneously,  the  total  amount  of  signal  information  may  be  significantly 
more  substantial  than  that  from  a  monostatic  radar.  With  the  signals  received  by  spatially 
separated  receiving  stations,  the  size,  form  and  relative  motion  of  a  target  can  be 
measured  with  higher  accuracy  and  in  a  shorter  time  interval.  Additionally,  two- 
dimensional  and  even  three-dimensional  radar  images  of  the  target  can  be  obtained  if  the 
transmitting  and  receiving  stations  are  spatially  coherent. 

5,  Increased  Resolution 

The  radar’s  resolution  capability  in  measuring  a  target’s  range,  bearing  or  velocity 
is  related  to  the  extent  of  the  response  of  echo  from  a  point  target.  High  resolution 
information  of  the  target  may  be  achieved  in  a  multistatic  radar  system  if  the  signals  at 
the  receiving  stations  are  spatially  correlated.  Using  the  hyperbolic  method,  these  signals 
undergo  correlation  processing  and  two  sources  of  mutually  uncorrelated  radiations, 
located  within  the  main  lobe  of  the  receiving  station’s  antenna,  can  be  resolved  by  the 
system. 


6,  Power  Advantages 

The  multitude  of  transmitting  and  receiving  stations  in  the  multistatic  radar 
system  has  additional  power  advantages  compared  to  monostatic  radar.  Each  receiving 
station  can  exploit  the  energy  transmission  from  all  transmitting  stations  and  enjoy 
significant  power  benefits.  When  the  baseline  distances  are  adequately  long,  scattered 
signal  fluctuations  are  statistically  independent  at  different  receiving  stations. 
Information  fusion  may  lead  to  additional  power  gain  due  to  fluctuation  smoothing. 
Additionally,  when  stations  are  separated  such  that  the  angle  between  directions  from  a 
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target  to  a  transmitting  and  a  receiving  station  approaches  180°,  the  radar  cross  section 
(RCS)  and  the  scattered  signal  intensity  may  increase  significantly  at  the  receiving 
stations.  This  increase  in  RCS  cannot  be  completely  reduced  by  stealth  technologies  such 
as  body  shaping  and  application  of  radar  absorbent  material  coating.  Furthermore,  since 
the  transmitter  and  receiver  are  spatially  separated,  receiver  protection  devices  such  as 
duplexer  and  circulator  are  not  required.  Consequently,  signal  power  losses  are  reduced. 

7.  Configurable  Coverage  Area 

System  geometry  and  fusion  algorithms  in  a  multistatic  radar  system  permit  the 
extension  of  coverage  area  in  the  required  direction.  If  the  radar  system  is  comprised  of 
mobile  transmitter  and/or  receiver  units,  the  reconfiguration  of  coverage  area  can  be 
achieved  easily. 

8,  Improved  Clutter  Rejection 

In  a  multistatic  radar  system,  since  the  transmitting  and  receiving  stations  are 
physically  separated,  the  intersection  volume  of  their  main  beams  may  be  much  less  than 
the  main  beam  volume  of  a  monostatic  radar.  Under  certain  conditions,  a  significant 
reduction  in  clutter  intensity  can  result.  Additionally,  moving  target  indication  (MTI) 
techniques  perform  better  in  a  multistatic  radar  than  a  monostatic  radar.  When  the  radial 
velocity  of  a  target  relative  to  a  monostatic  radar  is  zero,  MTI  techniques  are  useless.  In  a 
multistatic  radar  system,  a  moving  target  cannot  generally  present  zero  radial  velocity  to 
several  spatially  located  receiving  stations.  Similarly,  the  limitation  of  blind  radial 
velocities  in  a  monostatic  radar  is  overcome. 

I.  MAIN  DRAWBACKS 

With  the  increased  number  of  stations  and  components  in  a  multistatic  radar,  the 
added  complexities  also  impose  difficulties  in  the  implementation.  Some  of  these 
drawbacks  include  the  requirement  for  synchronization  and  mutual  alignment,  the  need 
for  direct  line  of  sight  between  stations  and  targets,  the  need  for  data  transmission  lines, 
and  increased  processing  requirements  [21]. 
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1.  Synchronization,  Phasing  and  Transmission  of  Reference  Signals 

Synchronization  between  the  transmitting,  reeeiving  stations  and  control  center  is 
necessary  for  accurate  target  coordinate  measurement  by  the  hyperbolic  methods.  When 
cooperative  signal  reeeption  is  used  in  the  multistatie  radar  system,  the  reeeivers  must 
know  the  signal  waveform  sent  out  by  the  transmitting  stations.  This  ean  be  aehieved  by 
signal  transmission  through  the  data  transmission  lines  or  transmission  of  synehronization 
eodes  to  eorreet  frequeney  and  signal  waveform.  For  eoherent  signal  proeessing,  a 
eommon  referenee  frequency  is  required  at  eaeh  station  to  couple  the  transmitters’  and 
receivers’  heterodyne  frequencies. 

2.  Requirement  for  Direct  Line  of  Sight  between  Stations  and  Targets 

The  coverage  of  a  multistatic  radar  system  is  limited  by  the  neeessity  for  direet 
lines  of  sight  between  stations  and  targets.  If  a  target  is  not  concurrently  visible  to 
several  transmitting  and  reeeiving  stations,  information  fusion  eannot  be  achieved.  This 
is  an  important  eonstraint,  espeeially  for  ground-based  multistatic  radar  systems  in  the 
detection  of  low-level  targets. 

3.  Need  for  Data  Transmission  Lines 

These  data  transmission  lines  are  used  for  signal  or  data  transmission  from  the 
stations  to  the  eontrol  eenter,  as  well  as  eommand  and  eontrol  information  transmissions. 
These  transmission  lines  need  to  be  protected  against  interferenee  and  direet  physieal 
attaeks. 


4,  Increased  Processing  Requirement 

In  a  multistatie  radar  system,  the  signifieant  inerease  in  information  from  a  target 
as  eompared  with  a  monostatie  radar  will  require  fast  proeessors  for  proeessing  these  data 
in  real  time.  Additionally,  the  eoordinate  eonversion  of  loeal  radar  data  from  eaeh 
reeeiving  station  into  a  eommon  eoordinate  system  and  the  interstation  data  assoeiation 
between  measurements  impose  added  proeessing  needs.  Furthermore,  geometrical  and 
tracking  algorithms  are  also  more  complex  than  those  applied  in  monostatic  radars. 
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5.  Need  for  Accurate  Station  Positioning  and  Mutual  Alignment 

The  accurate  fusion  of  target  information  from  several  spatially  located  stations 
requires  the  precise  knowledge  of  the  stations’  positions  and  their  alignment.  Positions  of 
stationary  stations  can  be  obtained  accurately  by  geodetic  methods.  However,  the 
determination  of  positions  of  mobile  stations  will  require  additional  navigation  methods 
such  as  GPS  information.  Errors  in  station  positions  and  their  misalignment  will  lead  to 
systematic  errors  of  target  location  determination. 

J,  IMPLEMENTATION  REQUIREMENTS 

It  is  important  that  signal  synchronization  and  spatial  coherence  among 
transmitting  and  receiving  stations  be  established. 

Signal  synchronization  is  required  between  the  transmitting  station  and  receiving 
station  for  range  measurement  [22].  Time  synchronization  can  be  accomplished  by 
receiving  the  signal  directly  from  the  transmitter.  This  transmitting  signal  can  be  sent 
directly  if  an  adequate  line  of  sight  exists  between  the  transmitter  and  receiver,  via  land 
line  or  other  communication  links  if  such  line  of  sight  cannot  be  achieved.  For  these 
direct  synchronization  means,  the  implementation  is  relatively  straightforward. 

Establishing  spatial  coherence  (amplitude  and  phase)  between  each  transmitting 
and  receiving  station,  and  the  coherent  combination  of  signals  from  all  receiving  stations 
are  essential  for  coherent  multistatic  radar  operation  [23].  For  a  ground-based  multistatic 
radar  system,  a  viable  solution  is  the  establishment  of  suitable  data  links  between  the 
elements  and  the  central  information  center. 

In  general,  a  multistatic  radar  could  be  easily  implemented  when  the  positions  of 
all  transmitting  and  receiving  stations  are  at  fixed  locations. 

For  the  purpose  of  this  thesis  research,  it  is  assumed  that  only  the  synchronization 
signal  and  Doppler  shift  information  received  by  each  receiving  station  are  available  to 
the  multistatic  radar  system  for  processing.  Thus,  the  challenge  is  to  resolve  the 
geometry  and  obtain  accurate  position  and  velocity  details  of  the  targets  based  on  limited 
amount  of  available  information. 
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V.  RADAR  MODEL  DESIGN 


A  baseline  design  and  approaeh  has  been  developed  to  demonstrate  the  eoneept  of 
high  resolution  imaging.  The  initial  radar  model  is  based  on  the  basie  multistatie  radar 
configuration  presented  in  the  preceding  chapter,  with  additional  parameters  included  to 
aid  in  the  derivation  of  the  required  parameters.  Figure  1 1  depicts  such  a  framework  and 
the  corresponding  parameters. 


T1 

Figure  1 1 .  Framework  for  initial  Multistatic  radar  model 

A.  ASSUMPTIONS 

The  radar  system  is  assumed  to  use  antennas  with  relatively  low  gain.  Fience,  the 
wide  beam  width  makes  it  unsuitable  for  direct  measurement  of  target  bearing.  It  is 
believed  that  the  implementation  of  such  a  system  is  feasible  and  the  focus  is  on  the 
localization  of  the  targets  and  the  resultant  formation  of  the  target  images. 


35 


For  the  purpose  of  this  thesis  researeh,  it  is  assumed  that  only  the  time 
synchronization  signal  and  Doppler  shift  information  reeeived  by  eaeh  reeeiving  station 
are  available  to  the  multistatie  radar  system  for  proeessing.  With  the  time 
synehronization  signal,  the  range  sum,  Rj  +  Rr,  eould  be  derived.  The  Doppler  shift  will 
result  in  the  knowledge  of  the  target’s  veloeity  veetor. 

Sinee  the  foeus  of  this  thesis  is  a  proof  of  coneept,  it  is  assumed  that,  for  ease  of 
demonstration,  the  operational  area  of  interest  is  a  flat  two-dimensional  plane.  After  the 
eoneept  has  been  proven,  the  more  realistie  three-dimensional  spaee  ean  be  extended  in  a 
relatively  straightforward  manner. 

The  target  is  assumed  to  behave  like  a  rigid  body.  Non-rigid  bodies  ean  be 
modeled  by  multiple  targets  with  varying  veloeity  vectors. 

A  flat  earth  has  also  been  assumed  to  simplify  the  derivation  and  demonstration. 

B,  MODEL  NOTATIONS 


The  notations  used  in  the  subsequent  derivation  are  presented  in  Table  2. 


Parameter 

Notation 

Description 

Bearing 

^Tgf^Tl-'^RO^Rl 

Heading  using  North-referenced  coordinate  system 

Bistatic  Angle 

Bistatic  angle  for  each  Transmitter-Target-Receiver 

bistatic  system 

Aspect  Angle 

^1,^2 

Aspect  angle  of  target  relative  to  bistatic  bisector  for 

each  Transmitter-Target-Receiver  bistatic  system 

Offset  Angle 

^off\  ’  ^offl 

Angle  that  baseline  for  each  Transmitter-Target- 

Receiver  bistatic  system  makes  with  horizontal  axis 

Subtended  Angle 

^R\  ’  ^R2  ’  ^T\R\ ’  ^T\R2 

Angle  subtended  within  each  bistatic  triangle 

Range 

RtI,  RrI,  Rr2 

Range  to  target 

Baseline  Distance 

Li,  L2 

Separation  between  Transmitter  and  Receiver  for  each 

Transmitter-Target-Receiver  bistatic  system 

Target  Velocity 

V 

Magnitude  of  target’s  velocity 

Table  2.  Notations  used  in  multistatic  radar  model 
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c. 


TARGET  BEARING  AND  RANGE  RELATIVE  TO  TRANSMITTER 


In  this  section,  the  objective  is  to  determine  the  position  of  the  target  relative  to 
the  transmitter.  Therefore,  the  parameters  to  be  determined  are;  target’s  bearing  relative 
to  the  transmitter  and  the  target’s  relative  range  to  the  transmitter,  Rji. 


From  Equation  (4.3),  the  derivation  oiRn  can  be  extended  in  the  multistatic  case 
to  include  one  for  each  of  the  two  bistatic  triangles. 


Rj 


{Rr+R^J-e 

2  )-  L  sin(^zi^  ) 


(4.3) 


For  the  multistatic  radar  model  (see  Figure  11),  the  baseline  for  each  transmitter-receiver 
pair  may  not  be  aligned  with  the  horizontal  axis;  hence,  additional  offset  angles  need  to 
be  considered  for  the  correct  derivation.  The  angle  (l)j  in  the  denominator  of  Equation 
(4.3)  will  need  to  include  and  6(,^2  • 


2  “  f-i  sin^^2,j  + 

_ {Rtx'^Rri)  _ 

2  +  R^2  ~  ^2 


(5.1) 


(5.2) 


With  the  availability  of  time  synchronization  signals,  A7),j  and  A7),2 ,  the  range 
sums,  Rti  +  Rri,  and  Rji  +  Rr2,  can  be  obtained  as 

Rjy  +  Rjii  —  cATj  +  Zyj  (5-3) 

+  Rr2  -  +  ^2  (5.4) 

Hence,  by  substituting  Equations  (5.3)  and  (5.4)  into  Equations  (5.1)  and  (5.2),  (l)^^  can 
be  estimated  by  the  numerical  solution  of 

_ (^n+^jji)  -L\ _ +i?2;2)  ~^2 _ _  Q 

2  “  f"!  sin^^2^j  +  2  +  ^r2  ^  ^2  +  0off2 
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Thereafter,  Rji  can  be  obtained  from  either  Equation  (5.1)  or  (5.2). 


D,  TARGET  BEARING  AND  RANGE  RELATIVE  TO  RECEIVERS 

Next,  the  positions  of  the  target  relative  to  the  receivers  need  to  be  determined. 
The  parameters  of  interest  are;  target’s  bearing  relative  to  the  receivers,  and 
the  target’s  relative  range  to  the  receivers,  Rri  and  Rr2.  After  Rji  and  ^zij-i  have  been 
determined,  Rri  and  Rr2  can  be  obtained  by  applying  the  cosine  rule  for  each  of  the 
bistatic  triangles.  Similar  to  the  derivations  for  the  transmitter,  the  offset  angles  due  to 
the  different  baseline  positions  have  to  be  considered. 


T1 

Figure  12.  Bistatic  triangle  involving  Tl,  R1  and  TGT 

For  the  bistatic  triangle  involving  Tl,  R1  and  TGT  (see  Figure  12),  the  derivation 
is  as  follows; 
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—  i?j.[  +  Zyj  —  2i?y,j_Zyj 


(5.6) 


Since 


,  therefore, 


i?^i  ■*■  2Rj.^L^  sin^^^j  +  0^^^  ^ 


Using  the  eosine  rule,  y^j  ean  be  determined  as 


y^i  =  eos  ' 


Rri  +  ~ 

2R,,Rt, 


Thereafter,  by  the  sine  rule,  6*^;  can  be  derived  as 


=  sin 


L, 


(5.7) 


(5.8) 


(5.9) 


Hence,  the  bearing  of  the  target  relative  to  receiver  R1  is 

^Rl=^Tl-^l  (5-10) 

It  should  be  noted  that  Equation  (5.10)  is  only  applicable  in  the  ease  when  the 
target  remains  on  the  same  side  of  the  extended  baseline  of  the  Transmitter-Reeeiver  pair. 
If  the  target’s  flight  path  is  sueh  that  it  crosses  to  the  other  side  of  the  extended  baseline, 
then  additional  correetions  will  need  to  be  done  to  obtain  the  eorreet  bearing.  This 
situation  of  target  traversing  to  the  other  side  of  the  extended  baseline  will  be  seen  when 
the  sum  of  the  following  angles  exceed  ti  /2  . 


When  (j)j^  +  +  dojji 


>  — ,  then 
2 


^R\  ^  ^T\R\  P\ 

^R\  ~  ^off  \  ~  ^R\ 

Similarly,  the  parameters  for  second  radar  receiver,  R2,  are  as  follows: 


^Ri  ~  yRri  7-2  2i?y,jZ,2  sin ^^2-1  +  ^offi ) 


(5.11) 
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1 


(5.12) 


=  cos 


p2  I  p2  _  t2 

^R2  ^T\  ^2 

2Rg^2^T\ 


^R2  ^T\  P2 


(5.13) 


E,  TARGET  VELOCITY  VECTOR 


The  target’s  velocity  information  can  be  extracted  in  two  ways.  Firstly,  it  can  be 
obtained  from  the  Doppler  shift  information  picked  up  by  the  receiving  stations. 
Secondly,  with  the  derived  target  position  information,  a  history  track  of  the  target  can  be 
maintained  and  the  corresponding  velocity  vector  can  be  derived.  The  required 
parameters  are  V,  (/)jg, ,  and  either  or  (5) .  The  Doppler  shift  method  will  be  discussed 

first.  The  amounts  of  Doppler  shift  observed  by  the  two  receiving  stations  are: 


(5.14) 

(5.15) 


By  geometry  (see  Figure  1 1), 


Substituting  into  Equation  (5.14)  and  equating  this  with  Equation  (5.15),  (5)  can  be 
found  by  the  least-squares  minimization  search  method  such  that  will  minimize  the 
following  function. 


/  = 


cos 


(5.17) 


hollowing  that,  the  magnitude  of  the  target’s  velocity  V  and  its  corresponding  velocity 
direction  are  derived. 
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(5.18) 


^Tgt  TC  (j)j^  82 


A 

2 


(5.19) 


It  should  be  noted  that  the  above  equations  might  produee  erroneous  results  when 
the  target  erosses  to  the  other  side  of  the  extended  baseline  of  the  Transmitter-Reeeiver 
pair.  This  ean  be  eorreeted  by  eheeking  the  magnitude  of  the  veloeity.  When  it  beeomes 
negative,  the  veloeity  veetor  direetion  will  also  have  n  radians  of  error. 
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VI.  RADAR  MODEL  IMPLEMENTATION  AND  RESULTS 


Based  on  the  model  introdueed  in  the  previous  ehapter,  a  proeedure  for 
implementation  will  be  deseribed  and  presented  in  this  ehapter.  The  seenario  is  a 
straightforward  test  ease  eomprising  one  transmitter  and  two  reeeivers.  The  results 
obtained  will  inelude  the  target’s  states  (loeation  and  veloeity)  at  eaeh  time  step.  These 
will  be  analyzed  and  the  errors  will  be  examined  in  detail. 

A,  SCENARIO  SETUP 

Consider  the  ease  of  a  simple  multistatie  radar  system  eomprising  just  one  radar 
transmitter  (Tl)  and  two  radar  reeeivers  (R1  &  R2),  with  two  airborne  targets  (Tgt  #1  & 
#2)  simulated  to  fly  in  different  flight  paths.  Tgt  #1  will  fly  in  a  sinusoidal  path  while 
Tgt  #2  will  fly  in  a  straight  line.  Figure  13  shows  a  sehematie  of  this  eonfiguration. 


;•  Tgt  #1 

/ 


Tgt  #2 


Figure  13.  Plan  view  sehematie  of  seenario  setup 

The  North  refereneed  eoordinate  system  will  be  adopted  for  the  seenario,  with  the 
origin  fixed  at  the  transmitter,  Tl.  The  eoordinate  positions  for  the  transmitting  and 
reeeiving  stations,  as  well  as  the  start  positions  of  the  airborne  targets  are  presented  in 
Table  3. 
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Entity 

Position 

Parameters 

Transmitter  (Tl) 

[0,  0] 

- 

Reeeiver  (Rl) 

[50,  20] 

Baseline  offset:  =  21.8° 

Reeeiver  (R2) 

[100,  0] 

Baseline  offset:  ~  0° 

Target  Initial  Position 

(Tgt#l) 

[0,  100] 

Veloeity:  Vix  =  100  m/s,  Viy  =100  m/s 

Target  Initial  Position 

(Tgt#2) 

[0,  100] 

Veloeity:  V2X  =100  m/s,  Viy  =  0 

Table  3.  Radar  stations  and  targets’  disposition  and  other  parameters 


In  this  setup,  different  target  flight  paths  and  different  radar  transmitter-reeeiver 
baselines  were  tested  for  eorreetness  of  information  extraetion.  T1  will  transmit  a  CW 
signal  at  a  fixed  frequeney  of  9  GHz  to  simulate  a  band-limited  system.  The  time  step 
used  is  1  ms  and  there  are  3000  steps  in  eaeh  simulation  run.  This  will  translate  to  3 
seeonds  of  target  observation  time. 

B,  TEST  DATA  GENERATION 

For  the  purpose  of  the  simulation,  test  data  is  divided  into  two  eategories:  data 
available  to  the  radar  system  and  data  for  results  eomparison  and  analysis. 

1.  Data  Available  to  Radar  System 

From  the  assumptions  for  this  thesis  researeh,  the  radar  system  will  only  have  two 
sets  of  data:  synchronization  time  delay  and  Doppler  shift.  In  the  simulation,  these  data 
will  be  generated  based  on  the  flight  paths  and  handed  over  to  the  radar  proeessing 
module  to  resolve  the  geometry  and  obtain  aeeurate  position  and  veloeity  details  of  the 
targets. 


2,  Data  for  Results  Comparison  and  Analysis 

The  data  for  results  eomparison  are  tabulated  in  Table  4.  These  data  are 
neeessary  to  model  eaeh  target  at  eaeh  time  step. 
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Data 

Description 

Target  Position 

Coordinates  of  the  target 

Target  Veloeity 

Magnitude  and  direetion  of  the  target 

Range  to  T1 

Range  from  target  to  transmitter  (Tl) 

Range  to  R1 

Range  from  target  to  reeeiver  (Rl) 

Range  to  R2 

Range  from  target  to  reeeiver  (R2) 

Bearing  from  T1 

North  refereneed  heading  of  target  from  transmitter  (Tl) 

Bearing  from  R1 

North  refereneed  heading  of  target  from  reeeiver  (Rl) 

Bearing  from  R2 

North  refereneed  heading  of  target  from  reeeiver  (R2) 

Bistatie  Angle  Tl-Rl 

Bistatie  angle  for  Tl-Target-Rl 

Bistatie  Angle  T1-R2 

Bistatie  angle  for  Tl-Target-R2 

Relative  Veloeity  from 

T1 

Veloeity  eomponent  along  bistatie  biseetor,  as  seen  by 

transmitter  (Tl) 

Relative  Veloeity  from 

R1 

Veloeity  eomponent  along  bistatie  biseetor,  as  seen  by 

reeeiver  (Rl) 

Relative  Veloeity  from 

R2 

Veloeity  eomponent  along  bistatie  biseetor,  as  seen  by 

reeeiver  (R2) 

Table  4.  Data  for  results  eomparison  and  analysis 


Based  on  the  test  ease,  a  Matlab  plot  of  the  flight  paths  and  the  radars  disposition 
are  as  shown  in  Figure  14. 
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Plan  View  of  Radar  Positions  and  Targets’  Flight  Paths 


Figure  14.  Radar  positions  and  target  flight  paths 


C.  LEAST  SQUARES  IMPLEMENTATION 

In  the  implementation  of  least  squares  analysis  during  radar  measurement,  the 
data  will  be  collected  in  blocks  of  5  time  steps  and  the  least  squares  curve  fit  will  be  done 
for  each  time  block.  Since  each  time  step  is  1  ms,  the  design  matrix  A  employed  initially 
to  solve  Ax  =  m  is  in  the  form 


1 

0.001 

0.001' 

0.001 

1 

0.002 

0.002' 

0.002 

1 

0.003 

0.003' 

0.003 

1 

0.004 

0.004' 

0.004 

1 

0.005 

0.005' 

0.005 
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While  it  may  appear  to  be  a  legitimate  design  matrix,  it  is  aetually  not  suitable  as 
it  involves  powers  of  t  when  it  is  a  small  number.  This  ean  be  verified  by  eheeking  the 
eondition  number  for  matrix 


cond(A)  =  5.893  x  10* 

It  is  elear  that  the  eondition  number  is  large  and  it  implies  numerieal  instability.  For 
example,  if  cond(A)  =  10^,  a  relative  data  error  of  the  order  of  10'^  may  imply  an  error  of 
100%  in  the  solution.  Henee,  matrix  A  is  badly  eonditioned  and  its  eolumns  are  nearly 
linearly  dependent.  This  behavior  is  eornmon,  and  it  illustrates  the  importanee  of  eareful 
measurement  aequisition.  A  better  basis  is  provided  by  powers  of  a  translated  and  sealed 
t  sueh  that 

(t- 0.003) 

0.002 


The  value  of  0.003  in  the  numerator  is  the  median  of  the  eoeffieients  used  for  the  eubie 
polynomial.  The  value  of  0.002  in  the  denominator  is  taken  to  be  the  differenee  between 
this  median  value  and  the  first  eoeffieient  value  (0.001).  The  aetual  design  matrix 
implemented  and  its  eondition  number  beeomes 


-1 

(-1)^ 

(-1)^ 

-0.5 

(-0.5)^ 

(-0.5> 

0 

0 

0 

0.5 

(0.5)^ 

(0.5)* 

1 

(1)^ 

(1)^ 

(6.1) 


cond(A)  =7.104 


The  resulting  design  is  well  eonditioned  and  it  is  used  for  the  analysis  of  position  and 
veloeity  information. 


D,  EXTRACTION  OF  TARGET  POSITION  INFORMATION 

The  analysis  of  the  signals  reeeived  by  the  radar  system  will  be  eategorized  into 
two  main  areas,  namely  derivation  of  position  information  from  the  synchronization  time 
delay  and  the  extraetion  of  veloeity  information  from  the  Doppler  shift.  The  derivation 

47 


of  target  position  information  will  involve  the  determination  of  transmitter  to  reeeiver 
“round-trip”  distanee,  target  bearing  and  range  relative  to  the  transmitter,  and  target 
bearing  and  range  relative  to  the  reeeivers. 

1.  Transmitter  to  Receiver  “Round-trip”  Distance 

From  the  synehronization  time  delays,  the  “round-trip”  range  from  transmitter  to 
target  to  reeeiver  was  first  derived  for  eaeh  time  step,  target  and  radar  transmitter-reeeiver 
pair.  The  range  was  obtained  by  multiplying  the  time  delay  with  the  speed  of  light  and 
adding  the  baseline  distanee  between  the  transmitter-reeeiver  pair. 

Range  Tl  Tgtl  Rl  =  Time  Delay  Tgtl  Rl  *  e  +  LI 

Range  Tl_Tgt2_Rl  =  Time  Delay  Tgt2_Rl  *  e  +  LI 

Range  Tl_Tgtl_R2  =  Time  Delay  Tgtl_R2  *  c  +  L2  ^  ^ 

Range  Tl_Tgt2_R2  =  Time  Delay  Tgt2_R2  *  e  +  L2 

For  the  test  ease,  the  two  targets  and  two  transmitter-reeeiver  pairs  will  involve  four  sets 
of  “round-trip”  ranges  for  time  step. 


2.  Target  Bearing  and  Range  Relative  to  Receivers  (R1  &  R2) 

With  the  target  bearing  and  range  relative  to  the  transmitter  known,  the  ranges  of 
the  target  relative  to  the  reeeivers  were  obtained  from  Equations  (5.7)  and  (5.11).  Next, 
the  bistatie  angles  were  derived  from  Equations  (5.8)  and  (5.12).  The  bearings  of  the 
target  relative  to  the  reeeivers  were  subsequently  obtained  from  Equations  (5.10)  and 
(5.13). 


=  COS  ' 


p2  I  p2  _  t2 

2RriRtI 


=  COS 


p2  I  p2  _  t2 

^R2  ^2 

2Rr2Rtx 


(5.7) 

(5.11) 

(5.8) 

(5.12) 
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Range  (North  Direction)  [m]  Range  (North  Direction)  [m] 


^R\  ^T\  A 


(5.10) 


^R2  ^T\  Pi  (5.13) 

Plots  of  the  measured  values  are  shown  in  Figures  15  and  16. 


Path  o1  Tgtl  as  seen  from  R1  (Measured) 


Path  of  Tgt2  as  seen  from  R1  (Measured) 


Range  (East  Direction)  [m] 


Figure  15.  Plot  of  Tgt  #1  and  Tgt  #2  from  measured  data  in  Reeeiver  R1 


Path  of  Tgtl  as  seen  from  R2  (Measured) 


Path  of  Tgt2  as  seen  from  R2  (Measured) 


'i 

5  100  - 


Range  (East  Direction)  [m] 


Figure  16.  Plot  of  Tgt  #1  and  Tgt  #2  from  measured  data  in  Reeeiver  R2 


In  view  of  the  laek  of  resolution  in  these  figures,  the  graphieal  results  appear  to  be 
aeeurate  enough  for  imaging  purposes.  Upon  detailed  eomparison  of  eaeh  data  point, 
measurement  errors  ean  be  ealeulated  and  analyzed.  The  average  errors  and 
corresponding  standard  deviations  of  the  measured  results  are  tabulated  in  Table  5. 
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Measurement  Errors  in 

Receiver  R1 

Receiver  R2 

Position 

Avg  Error 

Std  Dev 

Avg  Error 

Std  Dev 

Range  Error  for  Tgt  1  [m] 

-1.578x10-'^ 

1.727xl0'‘® 

-1.048x10'" 

2.214x10'" 

Range  Error  for  Tgt  2  [m] 

-4.597x10-'^ 

1.261x10'" 

1.445x10'" 

8.427x10'" 

Bearing  Error  for  Tgt  1  [rad] 

2.720x10'*^ 

2.060x10'" 

-7.915x10'" 

3.991x10'" 

Bearing  Error  for  Tgt  2  [rad] 

1.134x10'*^ 

3.261x10'" 

-2.722x10'" 

1.319x10'" 

Table  5.  Average  error  and  standard  deviation  in  target  position  errors 


E.  EXTRACTION  OF  TARGET  VELOCITY  INFORMATION 


The  extraetion  of  targets’  veloeity  information  will  involve  the  proeessing  of 
Doppler  information  established  by  reeeivers  R1  and  R2.  With  the  Doppler  information, 
the  least-square  minimization  seareh  method  was  adopted  in  Equation  (5.17)  to  obtain  the 
aspeet  angle  of  the  target  relative  to  the  bistatic  bisector. 


/  = 


cos 


(5.17) 


Thereafter,  the  magnitude  of  the  target’s  velocity  and  its  corresponding  velocity  direction 
were  derived  from  Equations  (5.18)  and  (5.19). 


V  = 


(5.18) 


^Tgt  —  ^  +  ^ri  ^2 


A 

2 


(5.19) 


Eigures  17  and  18  show  the  plots  of  the  velocity  magnitude  and  direction. 
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Direction  [radians]  Magnitude  [m/s] 


Figure  17.  Velocity  magnitude  of  Tgt  #1  and  Tgt  #2  from  measured  data 


Velocity  Direction  of  Tgtl  (Measured) 


Velocity  Direction  of  Tgt2  (Measured) 


50  100  150  200  250  300 

Range  (East  Direction)  [m] 


Figure  18.  Velocity  direction  of  Tgt  #1  and  Tgt  #2  from  measured  data 


From  these  figures,  it  may  appear  that  errors  are  more  pronounced  in  Tgt  #2  as 
compared  to  Tgt  #1.  In  reality,  the  errors  for  both  targets  are  comparable  and  the 
discrepancy  is  due  to  the  limitation  of  resolution  of  the  figures.  The  tabulated  errors  for 
velocity  magnitude  and  direction  are  tabulated  in  Table  6. 
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Measurement  Errors  in  Velocity 

Receiver  R1 

Avg  Error 

Std  Dev 

Magnitude  Error  for  Tgt  1  [m/s] 

-8.446X10''' 

0.492 

Magnitude  Error  for  Tgt  2  [m/s] 

-3.297x10'^ 

1.826 

Direction  Error  for  Tgt  1  [rad] 

5.834x10'^ 

3.174x10'^ 

Direction  Error  for  Tgt  2  [rad] 

5.836x10-^ 

3.190x10-^ 

Table  6.  Average  error  and  standard  deviation  in  target  veloeity  errors 
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VII.  ERROR  ANALYSIS 


A.  OVERVIEW 

The  errors  arising  as  a  result  of  the  measurements  ean  be  broadly  categorized  into 
four  areas,  namely  position  errors  in  range  and  bearing,  and  Doppler  errors  in  velocity 
magnitude  and  velocity  direction.  For  each  of  these  categories,  the  relevant 
measurements  are  compared  with  the  test  data  values,  and  the  differences  are  treated  as 
errors  and  plotted.  If  the  errors  are  significant,  they  will  be  manifested  in  the  radar  image 
as  imaging  artifacts.  Some  of  the  significant  errors  in  a  noiseless  environment  are  first 
analyzed.  They  serve  to  illustrate  the  effect  of  the  geometry  of  radar  systems  relative  to 
the  targets.  Although  the  measurements  can  be  considered  precise,  the  image  errors  can 
still  be  significant  if  the  condition  number  for  the  design  matrix  is  large.  Thereafter,  the 
errors  with  additive  Gaussian  noise  are  examined,  so  that  the  real  world  performance  of 
the  system  can  be  better  understood. 

B,  ERRORS  IN  BEARING  MEASUREMENT 

For  Receiver  Rl,  Figure  19  presents  the  errors  in  bearing.  It  is  evident  that  a 
“spike”  in  error  occurs  at  around  250  m.  This  occurrence  can  be  explained  by  reviewing 
the  flight  path  of  the  target  in  relation  to  the  disposition  of  the  radar  system.  The 
extended  baseline  of  the  transmitter-receiver  pair  for  Tl-Rl  is  as  shown  in  Figure  20.  It 
is  then  obvious  that  the  “spike”  occurs  at  the  intersection  of  this  extended  baseline  with 
the  flight  path  of  Tgt  #2.  This  is  reasonable  for  a  bistatic  radar  since  for  the  case  when 
the  target  lies  on  the  extended  baseline,  the  bistatic  angle,  =  0 ,  and  the  error  is 
maximum.  A  similar  “spike”  was  also  observed  in  the  errors  for  Tgt  #1. 
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Bearing  Error  of  Tgt2  as  seen  from  R1 
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Figure  19.  Bearing  error  for  Tgt  #2  as  seen  from  receiver  R1 


Extended  Baseline  of  T1  -R1 


Figure  20.  Extended  baseline  of  Transmitter-Receiver  pair  for  T 1  -R1 


54 


For  Receiver  R2,  Figure  21  presents  the  errors  in  bearing.  The  “spike”  occurs  for 
Tgt  #1  for  the  similar  reason  highlighted  for  Receiver  Rl.  Figure  22  depicts  the 
intersection  of  the  extended  baseline  for  Transmitter-Receiver  pair  T1-R2  with  Tgt  #Vs 
flight  path.  Such  “spikes”  are  not  evident  in  Receiver  R2  for  Tgt  #2  as  its  flight  path 
does  not  intersect  the  extended  baseline. 
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Bearing  error  for  Tgt  #1  as  seen  from  receiver  R2 


■  0 

Figure  21. 


X  10" 


Bearing  Error  of  Tgtl  as  seen  from  R2 
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Extended  Baseline  of  T1-R2 


Figure  22.  Extended  baseline  of  Transmitter-Receiver  pair  for  T1-R2 


C.  ERRORS  IN  VELOCITY  MEASUREMENT 

Figures  23  and  24  show  the  velocity  magnitude  errors.  For  Tgt  #1,  the  error 
“spikes”  occur  at  ranges  of  about  80  m  and  175  m.  The  first  “spike”  occurred  when  Tgt 
#1  was  turning  back  towards  the  South,  and  at  about  80  m  range,  its  aspect  angle  relative 
to  the  radar  system  was  90  degrees.  Consequently,  there  was  a  transition  of  the  radial 
velocity  through  zero  and  hence,  there  was  a  brief  moment  of  zero  Doppler  shift.  The 
second  “spike”  occurred  when  Tgt  #1  was  at  the  extended  baseline  of  the  Transmitter- 
Receiver  pair  Tl-Rl. 

For  Tgt  #2,  the  error  “spike”  happens  at  about  50  m.  This  is  a  special  case  where 
the  target  is  traveling  at  the  perpendicular  bisector  of  the  baseline  of  Transmitter- 
Receiver  pair  T1-R2  with  90  degrees  aspect  angle.  Consequently,  there  was  zero 
Doppler  shift  and  the  radar  system  responded  with  a  “spike”.  Figure  25  depicts  the 
erroneous  regions  for  velocity  magnitude  measurements. 
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Figure  23.  Velocity  magnitude  error  for  Tgt  #1 

Velocity  Magnitude  Error  of  Tgt2 
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Figure  24.  Velocity  magnitude  error  for  Tgt  #2 
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Erroneous  Areas  for  Velocity  Magnitude 


Figure  25.  Erroneous  regions  for  velocity  magnitude  measurements 


D.  ERRORS  WITH  ADDITIVE  GAUSSIAN  NOISE 

For  the  purpose  of  concept  demonstration,  the  errors  analyzed  so  far  have  been 
based  on  a  noiseless  environment.  In  actual  radar  operations,  noise  is  an  important  factor 
that  cannot  be  neglected.  Hence,  the  effect  of  noise  on  the  radar  model  was  tested  with 
additive  Gaussian  noise  applied  onto  the  synchronization  signal  and  the  Doppler  shifts 
measured  by  the  radar  system.  The  standard  deviation  used  was  two  percent  of  the 
measured  signal.  Figures  26,  27,  28,  and  29  present  samples  of  the  errors  for  Tgt  #1  with 
additive  Gaussian  noise. 
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Range  Error  of  Tgt1  as  seen  from  R1 
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Figure  26.  Range  error  for  Tgt  #1  as  seen  from  receiver  R1  (with  2%  Gaussian  noise) 
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Figure  27.  Bearing  error  for  Tgt  #1  as  seen  from  receiver  R1  (with  2%  Gaussian  noise) 
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28.  Velocity  magnitude  error  for  Tgt  #1  (with  2%  Gaussian  noise) 


Figure  29.  Velocity  direction  error  for  Tgt  #1  (with  2%  Gaussian  noise) 
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It  is  clear  that  the  errors  are  significantly  more  pronouneed  than  the  noiseless 
ease.  Furthermore,  sinee  the  order  of  magnitude  of  the  errors  are  larger,  the  error 
“spikes”  observed  in  the  noiseless  case  as  a  result  of  the  target  being  along  the  extended 
baseline  of  the  transmitter-receiver  pair  are  no  longer  apparent.  Nonetheless,  these 
results  are  still  accurate  enough  for  imaging  purposes. 

The  radar  model  was  tested  for  other  larger  values  of  additive  Gaussian  noise.  It 
was  noted  that,  when  the  additive  noise  was  inereased  to  ten  percent,  the  radar  model 
could  not  provide  much  conclusive  results.  Figures  30  and  31  show  the  sample  range  and 
bearing  errors.  This  was  due  to  the  use  of  the  functions  fzero  and  fminsearch  in  the 
Matlab  code.  In  the  simulation,  these  functions  encountered  extreme  values  and  eannot 
provide  reasonable  solutions  in  the  seareh  for  minimum  values.  The  algorithms  in  the 
Matlab  code  will  need  to  be  refined  to  handle  these  errors. 


Range  Error  of  Tgt1  as  seen  from  R1 
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Figure  30.  Range  error  for  Tgt  #1  as  seen  from  receiver  R1  (with  10%  Gaussian  noise) 
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Bearing  Error  of  Tgt1  as  seen  from  R1 
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Figure  3 1 .  Bearing  error  for  Tgt  #1  as  seen  from  reeeiver  R1  (with  10%  Gaussian  noise) 


Table  7  tabulates  the  standard  deviation  of  the  errors,  taken  to  be  the  root  mean 
square  deviation  from  the  mean,  as  a  result  of  various  amounts  of  additive  Gaussian 
noise. 


Gaussian 

Noise 

0.2% 

0.5% 

1.0% 

1.5% 

2.0% 

3.0% 

4.0% 

5.0% 

Range 

0.4122 

0.9378 

1.4966 

2.1441 

2.8388 

4.1039 

5.6526 

6.2225 

Bearing 

0.0242 

0.0551 

429.5414 

392.2657 

0.2024 

0.3073 

439.0524 

0.5210 

Veloeity 

Magnitude 

373.7993 

860.6418 

1.5707e3 

412.4362 

543.7135 

1.7291el3 

794.3593 

2.1248e4 

Veloeity 

Direetion 

0.3124 

0.4306 

429.5287 

392.2597 

9.5243+13 

1.8549el4 

3.5188el4 

3.3726el4 

Table  7.  Standard  deviation  of  the  different  types  of  errors  for  various  amounts  of 

Gaussian  Noise 
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Figures  32,  33,  34,  and  35  plot  these  errors  as  a  function  of  the  amount  of 
Gaussian  noise.  It  appears  that  the  relationship  between  the  standard  deviation  of  range 
errors  is  approximately  linear  to  the  amount  of  Gaussian  noise  added.  This  is  rational 
since  the  range  errors  are  largely  dependent  on  the  errors  in  time  synchronization  only. 
For  the  remaining  three  aspects,  namely  the  bearing,  velocity  magnitude  and  velocity 
direction  errors,  the  coupling  of  the  errors  in  the  time  synchronization  with  the  errors  in 
the  Doppler  frequencies  will  generally  result  in  larger  deviations  of  errors. 


Std  Dev  in  Range  Error  for  Various  Amounts  of  Gaussian  Noise 


Figure  32.  Standard  deviation  in  range  errors  for  various  amounts  of  Gaussian  noise 
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Std  Dev  in  Bearing  Error  tor  Various  Amounts  of  Gaussian  Noise 


Figure  33.  Standard  deviation  in  bearing  errors  for  various  amounts  of  Gaussian  noise 


std  Dev  in  Veiocity  Magnitude  Error  for  Various  Amounts  of  Gaussian  Noise 


Figure  34.  Standard  deviation  in  veloeity  magnitude  errors  for  various  amounts  of 

Gaussian  noise 
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Std  Dev  in  Velocity  Direction  Error  for  Various  Amounts  of  Gaussian  Noise 


Figure  35.  Standard  deviation  in  velocity  direction  errors  for  various  amounts  of  Gaussian 

noise 
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VIII.  RECOMMENDATIONS  AND  CONCLUSION 


A.  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  More  Complex  Scenarios 

The  multistatic  radar  system  analyzed  thus  far  was  based  on  the  straightforward 
disposition  of  one  transmitter  and  two  receivers.  Although  this  was  adequate  for  the 
proof  of  concept  of  providing  high-resolution  imagery  of  airborne  targets,  the 
employment  of  more  transmitting  and  receiving  stations  can  support  more  complex 
treatment  of  image  and  error  analysis.  With  more  transmitter-receiver  pairs,  additional 
algorithm  can  be  applied  to  eliminate  the  error  “spikes”  and  imaging  artifacts. 

Further,  the  scenario  was  tested  based  on  simple  two-dimensional  geometry.  A 
more  generalized  analysis  with  a  three-dimensional  perspective  could  reveal  issues  and 
challenges  that  degrade  the  performance  of  the  imaging  system.  For  instance,  due  to 
symmetry,  when  conducting  measurements  using  only  Doppler  information,  a  target 
trajectory  cannot  be  distinguished  from  its  mirror  image  about  a  vertical  plane  that  bisects 
the  baseline. 

Consequently,  as  the  scenario  becomes  more  complex,  the  algorithm  for  the  radar 
model  will  need  to  be  more  robust  to  handle  the  increased  complexity.  For  instance,  the 
limitations  in  the  use  of  Matlab  functions  /zero  and  fminsearch  need  to  be  overcome  to 
handle  scenarios  with  greater  noise  interference. 

2.  Performance  in  Real  World  Environment 

The  test  data  and  measurements  made  in  the  test  cases  were  simulated  to  be  in  a 
noise-free  environment.  The  inclusion  of  noise,  clutter,  and  other  environmental  effects 
such  as  multipath  in  the  simulations  will  allow  a  more  realistic  analysis  of  the  system’s 
performance.  Historically,  multipath  effects  have  prevented  the  exploitation  of  systems 
using  wide  antenna  beam  widths  for  precision  range  instrumentation  applications. 
Another  problem  unique  to  multistatic  sytems  is  multistatic  ghosts,  which  arise  when 
using  unaided  multilateration  techniques  without  transponders  in  a  multitarget 
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environment  [24],  These  false  target  indieations  are  generated  when  isorange  eontours 
from  two  targets  interseet  at  loeations  where  no  target  exists.  After  eonsidering  these  real 
world  effeets,  field  testing  of  sueh  a  system  will  eonstitute  an  important  eomponent 
towards  verifying  the  performanee  of  the  radar  system,  as  only  then,  ean  it  be  determined 
whether  this  concept  is  practicable. 

3,  Analysis  of  Bistatic  Ambiguity  Function 

Studies  have  shown  that  geometry  plays  an  important  role  in  the  shape  of  the 
bistatic  ambiguity  function  as  the  system  configuration  is  varied.  For  instance,  analysis 
in  [25]  has  proposed  algorithms  to  process  the  monostatic  ambiguity  function,  and 
thereafter,  extract  the  bistatic  factors.  When  the  bistatic  ambiguity  function  can  be 
determined  reasonably  well,  the  accuracy  of  the  multistatic  radar  system  can  be 
evaluated.  Additionally,  the  conditions  for  imaging  artifacts  can  be  better  understood, 
and  consequently,  the  artifacts  can  be  eliminated  more  effectively. 

B,  CONCLUSION 

The  multistatic  radar  offers  many  advantages  over  monostatic  radars  in  certain 
applications,  especially  since  the  receiving  stations  can  be  located  at  covert  and  distant 
sites  relative  to  the  transmitting  stations.  Furthermore,  CW  radars  are  relatively  simple 
and  inexpensive  to  employ  and  maintain.  Hence,  a  notional  CW  multistatic  radar  system 
for  high-resolution  imaging  was  examined.  This  concept  of  using  a  Doppler-only 
multistatic  radar  system  to  provide  high  resolution  imaging  of  airborne  targets  was 
demonstrated  to  be  attainable.  Through  an  understanding  of  conventional  imaging 
techniques  and  formulation  of  the  inverse  problem  in  radar  imaging,  the  demonstration 
radar  model  was  designed  to  acquire  an  accurate  position  and  velocity  of  the  targets.  The 
extraction  errors  resulted  from  the  range,  bearing  and  velocity  measurements  were 
congruent  with  the  physical  limitations  of  each  transmitter-receiver  pair.  Using  a 
multistatic  system,  the  geometrical  diversity  allowed  these  limitations  to  be  overcome. 
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APPENDIX 


%  Matlab  code  for  analysis  of 

%  High  Doppler  Resolution  Imaging  by  Multistatic  Continuous  Wave  Radars 
%  using  Constructive  Techniques 
"6 

%  by  Wei  Ting  Soh 
clear  all; 


f  =  9e9; 

"6 

frequency 

lambda  =  3e8  /  f; 

g, 

0 

wavelength 

dt  =  0.001; 

"6 

time  step 

Steps  =  3000; 

g, 

0 

number  of  steps 

%  Information  Available  to  Radar 


%  Radar  Positions 

Num_Tx  =1;  %  number  of  transmitters 

Num_Rx  =2;  %  number  of  receivers 

xy_Tx  =  [0  0 ] ; 
xy_Rx  =  [50  20;  100  0]; 

%  Angle  offset  for  diff  Tx_Rx  baseline  to  horizon, 

%  and  range  between  Tx  and  Rx 

theta_offset  =  zeros (Num_Rx) ;  %  offset  angle  for  each  baseline 

L  =  zeros (Num_Rx) ;  %  baseline  distance 

theta_offset  =  atan ( (xy_Rx ( : , 2 ) -xy_Tx (2 ) )  ./  (xy_Rx ( : , 1 ) -xy_Tx ( 1 ) ) ) ; 

L  =  sqrt ( (xy_Rx ( : , 2 ) -xy_Tx ( 1 , 2 ) ) . ^2  +  (xy_Rx ( : , 1 ) -xy_Tx ( 1 , 1 ) ) . ^2 ) ; 

%  Targets 

Num_Tgt  =2;  %  number  of  targets 

xyO  =  [0  100;  0  100]  ; 

V  =  [100  100;  100  0]  ; 


%  SET  UP  OF  TEST  CASE 


%  Data  Parameters 

data_xy_tgt  =  zeros (Steps, Num_Tgt , 2 ) ; 

each  time  step 

data_range_Tx  =  zeros (Steps, Num_Tgt, Num_Tx) ; 
at  each  time  step 

data_bearing_Tx  =  zeros (Steps, Num_Tgt, Num_Tx) ; 
Tx  at  each  time  step 

data_range_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 
at  each  time  step 

data_bearing_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 
Rx  at  each  time  step 

data_beta  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

targets  at  each  time  step 

data_velocity  =  zeros (Steps, Num_Tgt , 2 ) ; 
direction  of  targets  at  each  time  step 
data_Vr_Txl  =  zeros (Steps, Num_Tgt, Num_Tx) ; 

targets  to  Tx  at  each  time  step 


%  position  of  targets  at 
%  range  of  targets  to  Tx 
%  bearing  of  targets  to 
%  range  of  targets  to  Rx 
%  bearing  of  targets  to 
%  bistatic  angle  of 
%  velocity  magnitude  and 
%  relative  velocity  of 
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%  relative  velocity  of 


data_Vr_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

targets  to  Rx  at  each  time  step 

%  Information  available  to  radar  system 
time_delay  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

%  synchronization  time  delay  from  Tx  to  targets  to  Rx  at  each  time  step 

fd  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

%  Doppler  shift  of  targets  at  each  time  step 


for  k  =  l:Steps; 
t  =  k  *  dt; 

for  m  =  l:Num_Tgt; 

%  Target  position 

data_xy_tgt (k, m, : )  =  xyO (m, : )  +  v (m, : )  .*  [t, sin (2*t) ] ; 

%  Velocity 
if  k  >=  2 

[theta,  data_velocity (k, m, 1 ) ]  =  cart2pol ( (data_xy_tgt (k, m, 1 ) 
-  data_xy_tgt (k-1 , m, 1 ) ) /dt , (data_xy_tgt (k, m, 2 )  -  data_xy_tgt (k- 

l,m, 2) ) /dt) ; 


if  theta  <=  pi/2 

direction 


%  check  for  correct  sector 


data_velocity (k, m, 2 )  =  pi/2  -  theta; 
elseif  (theta  <=  1.5  *  pi)  and  (theta  >  pi/2) 
data_velocity (k, m, 2 )  =  theta  -  pi/2; 
elseif  (theta  <=  2  *  pi)  and  (theta  >  1.5*pi) 
data_velocity (k, m, 2 )  =  2.5*pi  -  theta; 

end; 

end; 


of  velocity 


for  n  =  l:Num_Rx; 


%  Relative  range 

data_range_Tx (k, m, 1 )  =  sqrt((data_xy_tgt(k,m,2)- 
xy_Tx ( 1 , 2 ) ) ^2  +  (data_xy_tgt (k, m, 1 ) -xy_Tx ( 1 , 1 ) ) ^2 ) ; 

data_range_Rx (k, m, n)  =  sqrt((data_xy_tgt(k,m,2)- 
xy_Rx (n, 2 ) ) ^2  +  (data_xy_tgt (k, m, 1 ) -xy_Rx (n, 1 ) ) ^2 ) ; 


%  Bearing  (measured  from  North) 

data_bearing_Tx (k, m, 1 )  =  atan ( (data_xy_tgt (k, m, 1 ) - 
xy_Tx (1,1) ) / (data_xy_tgt (k, m, 2 ) -xy_Tx (1,2) ) ) ; 


if  (data_xy_tgt (k, m, 2 ) -xy_Rx (n, 2 ) )  >=  0  %  Tgt  is  in  NE/NW 

quadrant 

data_bearing_Rx (k, m, n)  =  atan ( (data_xy_tgt (k, m, 1 ) - 
xy_Rx (n, 1 ) ) / (data_xy_tgt (k, m, 2 ) -xy_Rx (n,  2 ) ) ) ; 

elseif  (data_xy_tgt (k, m, 1 ) -xy_Rx (n, 1 ) )  >=  0  %  Tgt  is  in  SE 


quadrant 

data_bearing_Rx (k, m, n)  =  pi/2  + 
atan (abs ( (data_xy_tgt (k, m, 2 ) -xy_Rx (n, 2 ) ) / (data_xy_tgt (k, m, 1 ) - 
xy_Rx (n, 1) ) ) ) ; 

else  %  Tgt  is  in  SW  quadrant 

data_bearing_Rx (k, m, n)  =  -  pi/2  - 
atan (abs ( (data_xy_tgt (k, m, 2 ) -xy_Rx (n, 2) ) / (data_xy_tgt (k,  m,  1 ) - 
xy_Rx (n, 1) ) ) ) ; 

end; 


%  Beta  angle 
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data_beta (k, m, n)  =  acos ( (data_range_Tx (k, m, 1 ) ^2  + 
data_range_Rx (k, m, n) ^2  -  L(n)^2)  /  (2  *  data_range_Tx (k, m, 1 )  .* 

data_range_Rx (k, m, n) ) ) ; 

%  Time  delay  (between  Tx  and  Rx) 
time_delay (k, m, n)  =  (data_range_Tx (k, m, 1 )  + 

data_range_Rx (k, m, n)  -  L(n))  /  3.0e8; 

%  add  2%  Gaussian  noise  to  time  delay 
time_delay (k, m, n)  =  time_delay (k, m, n)  + 

0 . 02*time_delay (k,m, n) *randn; 

%  Relative  velocity  to  Tx 

data_Vr_Txl (k, m)  =  data_velocity (k, m, 1 )  *  cos  (pi  + 

data_bearing_Tx (k, m, 1 )  -  data_velocity (k, m, 2 ) ) ; 

%  Doppler  shift  for  Receivers 

data_Vr_Rx (k, m, n)  =  data_velocity (k, m, 1 )  *  cos (pi  + 

data_bearing_Tx (k, m, 1 )  -  data_velocity (k, m, 2 )  -  data_beta (k, m, n) ) ; 

fd(k,m,n)  =  (data_Vr_Txl (k, m)  +  data_Vr_Rx (k, m, n) )  /  lambda; 

%  add  2%  Gaussian  noise  to  Doppler  shift 
fd(k,m,n)  =  fd(k,m,n)  +  0 . 02*fd (k,m, n) *randn; 

end; 

end; 

end; 

%  Plot  scenario  setup 
%  figure; 

"0 

plot (xy_Tx ( 1 ) , xy_Tx (2 ) , ' kO ' , xy_Rx ( : , 1 ) , xy_Rx ( : , 2 ) , ' bX ' , data_xy_tgt ( : , 1 , 1 
) , data_xy_tgt (:,1,2), 'r', data_xy_tgt ( : , 2 , 1 ) , data_xy_tgt (:,2,2), 'b'); 

%  axis ( [-10  310  -20  210] ) ; 

%  xlabel (' Range  (East  Direction)  [m] ' ) ; 

%  ylabel (' Range  (North  Direction)  [m]  ' )  ; 

%  title  ('Plan  View  of  Radar  Positions  and  Targets''  Flight  Paths'); 

%  text (-3, -8,  ' T1  '  ) ; 

%  text  (47, 12,  'R1 ' ) ; 

%  text ( 97 , -8 , ' R2 ' ) ; 

%  text ( 70 , 1 90 , ' Tgt  #1'); 

%  text (75, 95, ' Tgt  #2'); 

%  Show  extended  baseline  for  Tl-Rl  pair 
%  figure; 

"6 

plot (xy_Tx ( 1 ) , xy_Tx (2 ) , ' kO ' , xy_Rx ( : , 1 ) , xy_Rx ( : , 2 ) , ' bX ' , data_xy_tgt ( : , 1 , 1 
) , data_xy_tgt (:,1,2), 'r', data_xy_tgt ( : , 2 , 1 ) , data_xy_tgt (:,2,2), 'b'); 

%  axis ( [-10  310  -20  210] ) ; 

%  xlabel (' Range  (East  Direction)  [m] ' ) ; 

%  ylabel (' Range  (North  Direction)  [m]  ' )  ; 

%  title (' Extended  Baseline  of  Tl-Rl'); 

%  text (-3, -8,  ' T1  '  ) ; 

%  text  (47, 12,  'R1 ' ) ; 

%  text ( 97 , -8 , ' R2 ' ) ; 

%  text ( 70 , 1 90 , ' Tgt  #1'); 

%  text (75, 95, ' Tgt  #2'); 

%  line([0  300],  [0  120 ],' Color ',' m ',' LineStyle ','-.') ; 

%  Show  extended  baseline  for  T1-R2  pair 
%  figure; 

"6 

plot (xy_Tx ( 1 ) , xy_Tx (2 ) , ' kO ' , xy_Rx {:,!), xy_Rx ( : , 2 ) , ' bX ' , data_xy_tgt ( : , 1 , 1 
) , data_xy_tgt (:,1,2), 'r', data_xy_tgt ( : , 2 , 1 ) , data_xy_tgt (:,2,2), 'b'); 
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%  axis ( [-10  310  -20  210] )  ; 

%  xlabel (' Range  (East  Direction)  [m] ' ) ; 

%  ylabel (' Range  (North  Direction)  [m] ' ) ; 

%  title (' Extended  Baseline  of  T1-R2  '  )  ; 

%  text (-3, -8,  ' T1  '  ) ; 

%  text  (47, 12,  'R1 ' ) ; 

%  text ( 97 , -8 , ' R2 ' ) ; 

%  text ( 70 , 1 90 , ' Tgt  #1'); 

%  text (75, 95, ' Tgt  #2'); 

%  line([0  300],  [0  0 ],' Color m LineStyle ; 

%  Show  extended  baseline  for  Tl-Rl  pair  &  perpendicular  bisector  for  Tl- 
R2  pair 
%  figure; 

"6 

plot (xy_Tx ( 1 ) , xy_Tx (2 ) , ' kO ' , xy_Rx ( : , 1 ) , xy_Rx ( : , 2 ) , ' bX ' , data_xy_tgt ( : , 1 , 1 
) , data_xy_tgt (:,1,2), 'r', data_xy_tgt ( : , 2 , 1 ) , data_xy_tgt (:,2,2), 'b'); 

%  axis ( [-10  310  -20  210] ) ; 

%  xlabel (' Range  (East  Direction)  [m] ' ) ; 

%  ylabel (' Range  (North  Direction)  [m] ' ) ; 

%  title  (' Erroneous  Areas  for  Velocity  Magnitude'); 

%  text (-3, -8,  ' T1  '  ) ; 

%  text  (47, 12,  'R1 ' ) ; 

%  text ( 97 , -8 , ' R2 ' ) ; 

%  text ( 70 , 1 90 , ' Tgt  #1'); 

%  text (75, 95, ' Tgt  #2'); 

%  line([0  200],  [0  80 ],' Color m LineStyle ; 

%  line([50  50],  [-10  120 ],' Color m LineStyle ; 

%  line([78  78],  [170  2 10 ],' Color m LineStyle ; 

"6 

%  DATA  PROCESSING  AND  SOLUTIONING 

x_scale  =  [0.1:0.1:300] ';  %  set  correct  scale  for  x-axis  when  plotting 

range_Tx  =  zeros (Steps, Num_Tgt, Num_Tx) ; 

at  each  time  step 

bearing_Tx  =  zeros (Steps, Num_Tgt, Num_Tx) ; 

Tx  at  each  time  step 

range_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

at  each  time  step 

bearing_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

Rx  at  each  time  step 

range_Tx_Rx  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

Tx  to  Tgt  to  Rx  at  each  time  step 

"0 

%  round  trip  range  from  Tx  to  Tgt  to  Rx 

for  m  =  l:Num_Tgt; 

for  n  =  l:Num_Rx; 

%  range  of  Tl-Tgt-Rl,  Tl-Tgt-R2 

range_Tx_Rx ( : , m, n)  =  time_delay ( : , m, n)  *  3.0e8  +  L(n); 


%  range  of  targets  to  Tx 
%  bearing  of  targets  to 
%  range  of  targets  to  Rx 
%  bearing  of  targets  to 
%  round  trip  range  from 


end; 

end; 

"0 

%  for  Tx 

%  output:  range_Tx (k, m, n)  &  bearing_Tx (k, m, n) 
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for  m  =  l:Num_Tgt; 
for  k  =  l:Steps; 

thetaT  =  @ (x)  (range_Tx_Rx (k, m, 1 ) ^2-L ( 1 ) ^2 )  /  (2  * 

(range_Tx_Rx (k, m, 1 )  -  L(l)*sin(x  +  theta_of f set ( 1 ) ) ) )  - 

(range_Tx_Rx (k, m, 2 ) ^2-L (2 ) ^2 )  /  (2  *  (range_Tx_Rx (k, m, 2 )  -  L(2)*sin(x  + 

theta_offset (2) ) ) )  ; 

bearing_Tx (k, m, 1 )  =  f zero (thetaT,  0 )  ; 

range_Tx (k, m, 1 )  =  (range_Tx_Rx (k, m, 1 ) ^2-L  ( 1 ) ^2 )  /  (2  * 

(range_Tx_Rx (k, m, 1 )  -  L ( 1 ) *sin (bearing_Tx (k, m, 1 )  +  theta_of f set ( 1 ) ) ) )  ; 

end; 

end; 


%  for  all  Rx 

%  output  -  position:  range_Rx (k, m, n)  &  bearing_Rx (k, m, n) 

%  velocity:  velocity (k, m, n, 2 )  [magnitude  &  bearing] 

beta  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

targets  for  each  Tx-Rx  pair  at  each  time  step 
delta  =  zeros (Steps, Num_Tgt, Num_Rx) ; 

for  each  Tx-Rx  pair  at  each  time  step 
velocity_dop  =  zeros (Steps, Num_Tgt , 2 ) ; 
direction  of  targets  at  each  time  step 

for  m  =  l:Num_Tgt; 

for  n  =  l:Num_Rx; 

for  k  =  l:Steps; 

range_Rx (k, m, n)  =  sqrt (range_Tx (k, m, 1 ) ^2  +  L(n)^2  -  2  * 
range_Tx (k, m, 1 )  *  L(n)  *  sin  (bearing_Tx (k, m, 1 )  +  theta_offset (n) ) ) ; 

beta(k,m,n)  =  acos ( (range_Rx (k, m, n) ^2  +  range_Tx (k, m, 1 ) ^2  - 
L(n)^2)  /  (2  *  range_Rx (k, m, n)  *  range_Tx (k, m, 1 ) ) ) ; 

%  angle  Tgt-Tx_Rx 

theta_Tx_Rx  =  asin (range_Rx (k, m, n)  *  sin (beta (k, m, n) )  / 

L(n)  )  ; 

if  bearing_Tx (k, m, 1 )  +  theta_Tx_Rx  +  theta_offset (n)  >  pi/2 

+  0.000001 

theta_Rx  =  pi  -  theta_Tx_Rx  -  beta(k,m,n); 

bearing_Rx (k, m, n)  =  1.5*pi  -  theta_offset (n)  -  theta_Rx; 

else 

bearing_Rx (k, m, n)  =  bearing_Tx (k, m, 1 )  -  beta(k,m,n); 

end; 

%  velocity  relative  to  each  receiver 

d2  =  @ (x)  (fd(k,m,l)  *  cos (x)  *  cos (beta (k, m, 2 ) /2 )  - 

fd(k,m,2)  *  cos ( 0 . 5* (beta (k, m, 2 ) -beta (k, m, 1 ) ) +x)  * 

cos (beta (k,m, 1) /2) ) ^2; 

delta(k,m,2)  =  fminsearch (d2 ,  1 )  ; 

delta(k,m,l)  =  delta(k,m,2)  +  (beta (k, m, 2 ) -beta (k, m, 1 ) ) /2 ; 

velocity_dop (k, m, 1 )  =  fd(k,m,2)  *  lambda  /  (2  * 

cos (delta (k, m, 2 ) )  *  cos (beta (k, m, 2 ) /2 ) ) ; 

velocity_dop (k, m, 2 )  =  pi  +  bearing_Tx (k, m, 1 )  -  delta(k,m,2) 

-  beta (k, m, 2 ) /2 ; 


%  Bistatic  angle  of 
%  Delta  angle  of  targets 
%  velocity  magnitude  and 
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%  process  negative  velocities 
if  velocity_dop (k, m, 1 )  <  0 

velocity_dop (k, m, 1 )  =  abs (velocity_dop (k, m, 1 ) ) ; 
if  velocity_dop (k, m, 2 )  >  pi 

velocity_dop (k, m, 2 )  =  velocity_dop (k, m, 2 )  -  pi; 

end; 

end; 

end; 

end; 

end; 


%  Least  square  processing:  Ax  =  B 
%  Error  calculation 


window  =  5; 

soln_Range  =  []; 

soln_Bearing  =  []; 

Processed_Range  =  []; 
Processed_Bearing  =  []; 
Processed_VelMag  =  []; 
Processed_VelDir  =  []; 

error_Range  =  [  ]  ; 

error_Bearing  =  []; 

error_VelMag  =  []; 

error_VelDir  =  []; 


%  Design  matrix 
A  =  []  ; 

for  s  =  1: window; 

tt  =  (s*dt  -  0.003)  /  0.002; 
design  matrix 

A  =  [A;  1  tt  tt^2  tt^3] ; 

end; 


%  time  scaling  for  better  basis  of 
%  design  matrix 


for  m  =  l:Num_Tgt; 

for  n  =  l:Num_Rx; 

soln_Range  =  []; 
soln_Bearing  =  []; 
soln_VelMag  =  []; 
soln_VelDir  =  []; 

for  1=1: (window) :3000;  %  divide  entire  range  into  blocks  of 

' window ' 

B1  =  [ ] ;  %  A  *  xl  =  B1 
B2  =  [  ]  ;  %  A  *  x2  =  B2 
B3  =  []  ; 

B4  =  []  ; 

for  s  =  1: window;  %  extract  data  for  the  window 


B1 

= 

[Bl;  range_Rx (l+s-l,m, n) ] ; 

B2 

[B2;  bearing_Rx (l+s-l,m, n) ] ; 

B3 

[B3;  velocity_dop ( 1+s-l , m, 1 ) ] ; 

end 

B4 

f 

[B4;  velocity_dop ( 1+s-l , m, 2 ) ] ; 

xl 

=  A 

★ 

Isqr (A,  Bl,  le-16); 

x2 

=  A 

★ 

Isqr (A,  B2,  le-16); 

x3 

=  A 

•k 

isqr (A,  B3,  le-16); 

x4 

=  A 

•k 

isqr (A,  B4,  le-16); 
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soln_Range  =  [soln_Range;  xl]; 
soln_Bearing  =  [ soln_Bearing;  x2]; 
soln_VelMag  =  [ soln_VelMag;  x3]  ; 
soln_VelDir  =  [ soln_VelDir ;  x4]; 

end; 

Processed_Range ( : , m, n)  =  soln_Range; 

Processed_Bearing ( : , m, n)  =  soln_Bearing; 

error_Range ( : , m, n)  =  soln_Range  -  data_range_Rx ( : , m, n) ; 
error_Bearing ( : , m, n)  =  soln_Bearing  -  data_bearing_Rx ( : , m, n) 

%  Receiver  range/bearing  errors 
figure; 

plot (x_scale, error_Range ( : , m, n) ) ; 

title ([' Range  Error  of  Tgt ' , int2str (m) ,  '  as  seen  from 

R ' , int2str (n) ] ) ; 

xlabel (' Range  (East  Direction)  [m] ' ) ; 
ylabel ( 'Absolute  Error  [m]  ' )  ; 

figure; 

plot (x_scale, error_Bearing ( :  ,  m,  n) ) ; 

title ([' Bearing  Error  of  Tgt ' , int2str (m) , '  as  seen  from 
R ' , int2str (n) ] ) ; 

xlabel (' Range  (East  Direction)  [m] ' ) ; 
ylabel ( 'Absolute  Error  [rad]  ' )  ; 


end; 

Processed_VelMag ( : , m)  =  soln_VelMag; 

Processed_VelDir ( : , m)  =  soln_VelDir; 

error_VelMag ( : , m)  =  soln_VelMag  -  data_velocity ( : , m, 1 ) ; 

error_VelDir ( : , m)  =  soln_VelDir  -  data_velocity ( : , m, 2 ) 


%  Doppler  errors 
figure; 

plot (x_scale, error_VelMag ( : , m) ) ; 

title ([ 'Velocity  Magnitude  Error  of  Tgt ' , int2str (m) ] ) ; 
xlabel (' Range  (East  Direction)  [m] ' ) ; 
ylabel (' Absolute  Error  [m/s]'); 

figure; 

plot (x_scale, error_VelDir ( : , m) ) ; 

title ([ 'Velocity  Direction  Error  of  Tgt ' , int2str (m) ] ) ; 
xlabel (' Range  (East  Direction)  [m] ' ) ; 
ylabel ( 'Absolute  Error  [rad]  '  )  ; 


end; 
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